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Summary. — 

After the discovery that supermassive black holes (SMBHs) are ubiquitous at the 
center of stellar spheroids and that their mass Mbh, in the range 
is tightly related to global properties of the host stellar system, the idea of the 
co-evolution of elliptical galaxies and of their SMBHs has become a central topic 
of modern astrophysics. Here, I summarize some consequences that can be derived 
from the galaxy scaling laws and present a coherent scenario for the formation and 
evolution of elliptical galaxies and their central SMBHs, focusing in particular on 
the establishment and maintenance of their scaling laws. In particular, after a first 
observationally based part, the discussion focuses on the physical interpretation of 
the Fundamental Plane. Then, two important processes in principle able to destroy 
the galaxy and SMBH scaling laws, namely galaxy merging and cooling flows, are 
analyzed. Arguments supporting the necessity to clearly distinguish between the 
origin and maintenance of the different Saling Laws, and the unavoidable occurrence 
of SMBH feedback on the galaxy Interstellar Medium in the late stages of galaxy 
evolution (when elliptical galaxies are sometimes considered as "dead, red objects"), 
are then presented. At the end of the paper I will discuss some implications of the 
recent discovery of super-dense ellipticals in the distant Universe. In particular, I 
will argue that, if confirmed, these new observations would lead to the conclusion 
that at early epochs a relation between the stellar mass of the galaxy and the 
mass of the central SMBH should hold, consistent with the present day Magorrian 
relation, while the proportionality coefficient between Mbh and the scale of velocity 
dispersion of the hosting spheroids should be significantly smaller than that at the 
present epoch. 
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1. Introduction 

The mutual interplay between supermassive black holes (hereafter SMBHs) at the 
center of stellar spheroids ( 1 ) and their host systems is now established beyond any rea- 
sonable doubt, as indicated by the remarkable correlations found between host galaxy 
properties and the masses of their SMBHs (e.g. [3]- [10]). More specifically, it is now 
believed that all early-type galaxies with Mq< — 18 mag ([11]) host a central SMBH 
(e.g. [12]- [14]), whose mass Mbh scales linearly with the spheroid stellar mass M*; the 
correlation of Mbh with the central stellar velocity dispersion a Q of the host galaxy is 
even tighter. It is then natural to argue (e.g. [15]- [27]) that the central SMBHs have 
played an important role in the processes of galaxy formation and evolution, the imprint 
of which is represented by the Scaling Laws (hereafter SLs) mentioned above. As an addi- 
tional supporting argument, several groups have noted the link between the cosmological 
evolution of QSOs and the formation history of galaxies (e.g. [28]- [32], see also [33]- [37]). 

In addition to the SLs of their central SMBHs, early-type galaxies are also known to 
follow well defined empirical SLs relating their global observational properties, such as 
total luminosity L, effective radius R e , and central velocity dispersion a . Among others 
we recall the Faber & Jackson ([38], hereafter FJ), the Kormendy ([39]), the Fundamental 
Plane ([40, 41], hereafter FP), the color-cr ([42]), and the Mg 2 -er (e.g. [44]-[46]) relations. 

Clearly, all together these scaling relations reveal the remarkable homogeneity of early- 
type galaxies, provide invaluable information about their formation and evolution, and 
set stringent requirements that must be taken into account by any proposed galaxy- 
SMBH formation scenario. We recall that two alternative scenarios for the formation of 
Es have been proposed. In the monolithic collapse picture, ellipticals are formed at early 
times by dissipative processes (e.g., [47]- [49]; see also [50, 51]), while in the hierarchical 
merging scenario spheroidal systems are the end-products of several merging processes of 
smaller galaxies, the last major merger taking place in relatively recent times, i.e. at z<l 
(e.g. [52]- [56]). Each of the two scenarios scores observational and theoretical successes 
and drawbacks (e.g. [57, 58]). For example, the merging picture (in its dry flavor, i.e. 
neglecting the role of the dissipative gas), could be supported by some observational 
data suggesting that a fraction of red galaxies in clusters at intermediate redshift are 
undergoing merging processes: these galaxies could be the progenitors of present-day 
early- type galaxies (e.g. [59]). At the same time, it is not clear how repeated merging 
events can produce a class of objects following striking scaling laws involving their global 
structure, dynamics, and stellar population properties, while theoretical investigations 
showed that the dynamical processes expected to follow the strongly dissipative phases 
of monolithic collapse apparently lead to systems surprisingly similar to real Es (e.g., 
see [60, 61]). 

In this paper I will summarize some selected topics of research in this field and present 
a possible coherent scenario for the co-evolution of SMBHs and their host spheroids. 
Because of the enormous body of dedicated literature now available, I will focus my 
attention on a few selected topics. Very important arguments, such as the dynamics 
of binary SMBHs and their effects on the central regions of galaxies (e.g. [62]- [65] and 



( 1 ) The term early-type galaxies is generically used for galaxies belonging to the family of 
elliptical galaxies (Es), SO galaxies, dwarf ellipticals (dE), and dwarf spheroidals (dSph); the 
class of stellar spheroids is made of early-type galaxies and bulges of spiral galaxies. For a 
detailed account of the observational properties of these classes see, e.g., [1, 2]. 



CO-EVOLUTION OF ELLIPTICAL GALAXIES AND THEIR. CENTRAL BLACK HOLES 



3 



references therein), or the growth of SMBHs via accretion of stars (e.g. [66] and references 
therein) are just mentioned here but not discussed. The paper is organized as follows. 
After a first, observationally based part (Sect. 2), where the main SLs followed by early- 
type galaxies and by their central SMBHs are described, in Sect. 3 I will focus on the 
physical interpretation of the FP: this phenomenological part is a prerequisite to the 
construction of a possible formation scenario. Then (Sect. 4) I discuss two important 
physical mechanisms in principle able to destroy the galaxy and SMBH SLs, namely 
galaxy merging and cooling flows. In this context, I present some arguments that require 
that we distinguish between the origin and maintenance of the SLs and the unavoidable 
occurrence of SMBH feedback on the galaxy Interstellar Medium in the late stages (when 
Es are sometimes considered as "dead, red objects"). This last point is strictly related 
to the solution of the long-standing problem of cooling flows in Es (and clusters), and 
to the quiescence of Active Galactic Nuclei (AGN) in low-redshift systems. Section 5 
addresses the problem of the origin of the SLs, while Sect. 6 provides a summary of the 
main results, with a short discussion of some very recent observational findings, i.e. the 
surprisingly high stellar density of Es at high redshift. 



2. Basic observational facts 

In this Section I review some of the most important SLs of early-type galaxies 
(Sect. 2.1), of their central SMBHs (Sect. 2.2), and the counterpart of the FP, FJ, and 
Kormendy relations followed by galaxy clusters (Sect. 2.3). Finally, the AGN underlu- 
minosity problem is presented (Sect. 2.4). 

2T. Scaling Laws of early-type galaxies . - 

2T.1. The Fundamental Plane. Early-type galaxies can be characterized by three 
main observable global scales: the circularized effective radius R e = \fa e b e (where a e 
and b c are the major and minor semi-axis of the effective isophotal ellipse, i.e., the 
ellipse containing half of the projected system luminosity), the central projected velocity 
dispersion er (often referred to an aperture radius of R e /8, e.g. [67]), and the mean 
effective surface brightness within R e , SB e = — 2.51og(/) c (where (/) c = Le/27ri? 2 , and 
Lb is the luminosity of the galaxy, for example in the Johnson B-band). It is well 
known that Es (with a — 100 km s -1 ) do not populate uniformly this three-dimensional 
parameter space; rather, they are confined to the vicinity of a narrow logarithmic plane 
([40, 41]) 

(1) logi? c = a log do +I3SB c + j 

thus called the Fundamental Plane (FP). The coefficients a, (3, and 7 depend slightly on 
the photometric band considered. By measuring R c in kpc, o~ a in km s , and SB C = 
42.0521 — 2.51og(L B /27ri? 2 ) in mag/arcsec 2 , where L B is expressed in units of the solar 
Blue luminosity, reported values are a = 1.25 ± 0.1, (3 — 0.32 ± 0.03, 7 = —8.895 
(e.g [68]-[74]). An alternative expression for the FP has been obtained by [75], with 
the introduction of the k coordinate system, in which the new variables are a linear 
combination of the observables: 

(2) fcl = ___ k2= _ , k3= _ . 
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Fig. 1. - The distribution of Virgo (closed boxes) and Coma (crosses) ellipticals in the (fci, k,2, £3) 
space. The upper panel shows the FP edge-on; in the lower panel the FP is seen nearly face-on 
(from [94]). 



In particular, when projected on the (fci, £3) plane, the FP is seen almost edge-on and it 
is considerably thin, while the distribution of galaxies in the (k\, fe) plane is considerably 
broader. For example, Virgo ellipticals are distributed on the (fci, £3) plane according to 
the best-fit relation 

(3) k 3 = 0.15fci +0.36 

(when adopting respectively, kpc, kms -1 and L-qq pc -2 as length, velocity and surface 
brightness units, based on a Virgo distance of 20.7 Mpc, see Fig. [1]). The systematic 
increase of kz along the FP described by eq. (3) and the nearly constant and very small 
dispersion of fc 3 at every location on the FP (for example with <r[fc 3 ] ~ 0.05 for Virgo 
ellipticals) are usually referred to as tilt and thickness of the FP. 

21.2. The Faber- Jackson and the Kormendy relations. Connected with the FP are the 
less tight FJ and Kormendy relations. These scaling relations have often been considered 
just as projections of the FP with no additional information, and, for this reason, their 
constraining power on galaxy formation scenarios has been underestimated. However, 
this is not fully correct, because these two relations describe where Es are distributed on 
the FP and so, even if characterized by a larger scatter than the edge-on view of the FP, 
they contain important information on the galactic properties. The proposed original 



CO-EVOLUTION OF ELLIPTICAL GALAXIES AND THEIR. CENTRAL BLACK HOLES 



5 



form of the FJ relation was Lb oc it™, with n w 4, while [76] found that the double-slope 
fit 

(4) ^ * 0.23 ( Oo A " + 0.62 f g ° r ) 4 ' 2 
V ; 10 Ibg VSOOkms- 1 / VSOOkms- 1 / 

provides a better description of the data. Note that for small galaxies and bulges 
(<t ^ 170 km s -1 ) the single power-law fit would give n ~ 2.4, considerably smaller than 
4 (see also [77, 78] for the case of the galaxies of the Virgo and Coma clusters). Re- 
cent measurements, based on the large data set obtained from the Sloan Survey SSDS 
([79]), converge to an exponent 4 in the K band; thus, this value is currently adopted in 
applications of the FJ relation to high-luminosity galaxies. 

The total luminosity of bright spheroidal systems also correlates with their length 
scale as measured by R e : in fact, such Kormendy relation can be written in the form 

(5) R c oc L B a , 

where the exponent a is strongly dependent on the galaxy sample used, and is found in 
the range 0.88 <a< 1.62 (e.g. [80]). The latest estimates appear to converge to a value 
a ~ 0.7 or less, as a function of waveband ([79]). 

2T.3. Structural weak homology. The empirical R 1/>4 luminosity "law" ([81], see eq. [6] 
below with m = 4), has long been recognized to fit the surface brightness profiles I(R) 
of Es successfully, to the point that Es are routinely identified by means of this char- 
acteristic photometric signature. It has no free parameters and depends on two well 
defined physical scales, the effective radius R c and the central surface brightness Jo- 
In practice the overall i? 1 / 4 fit is characterized by residuals typically of the order of 
0.1 — 0.2 mag/arcsec 2 (e.g. [82]-[86]). These deviations from the i? 1 / 4 law, although small, 
are often larger than the typical observational errors involved. In fact, the surface bright- 
ness distribution of ellipticals is better described by the Sersic ([87]) i? 1 / m -law 



(6) I(R) = I exp 



-Km) ( x 



R x 



lo 



R 2 C 27rmr(2m) 



(e.g. [88]-[100]). In the above equation T is the complete Gamma function. An extremely 
accurate analytical representation of the factor b(m) for ra > 1 is given by its truncated 
asymptotic expansion b{m) ~ 2m — 1/3 + 4/405™ ([101]. For a dynamical analysis of 
models with R x / m projected density profile see, e.g. [102]-[105]; see also [106]). 

It has soon become clear that the shape parameter m of the R}l m law correlates with 
global quantities such as total luminosity and effective radius (e.g. [90, 97], [107]- [111]), 
with m increasing with luminosity from ~ 1/2 up to ~ 15. For example [112, 113] report 
the following relations 

(7) logm- 0.28 + 0.52 logi? c ; m ~ -19.082 + 3.0275 logi. 

This remarkable luminosity dependence of the surface brightness profile has been called 
weak homology. 

In addition to the global trends described above, also local relations - relevant in the 
present context - have been found. For example, ground based observations ([114]) and 
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Fig. 2. - R 1 ^" 1 fit. R is in arcsec, fi in mag/arcsec 2 . The data points are in the Blue band, taken 
from [112, 113]. Note the limiting surface brightness reached for each galaxy (from [100]). 



Hubble Space Telescope data show that the volume luminosity profiles of Es approaches 
the power-law form p(r) oc r~ 7 at small radii, with < 7 < 2.5 ([115]-[121]). In 
particular, HST observations reveal that at small radii some profiles are rather flat (core 
galaxies) while others are characterized by steep cusps (power-law galaxies) . In general, 
core profiles are common among bright Es, but fainter systems tend to have power- 
law cusps; remarkably, other galaxy global properties are related to the presence of the 
core ([122, 123]). In the context of surface brightness profiles, it has also been found 
that while the surface brightness profiles of power-law galaxies are well fitted by the 
R}l m law all the way into the center, within a break radius Rb the surface brightness 
profile of core galaxies stays well below their global best-fit Sersic model; in addition, 
intermediate and low-luminosity galaxies often present additional star clusters at their 
center (e.g. [Ill, 124]-[128]). From these central-to-global relations, very interesting 
consequences can be derived (see Sects. 5 and 6). 

2 "2. Scaling laws of the central SMBHs. - The masses Mbh of central SMBHs lie 
in the range 10 6 — 10 9 M Q and correlate surprisingly well with several global and local 
properties of the host systems (e.g. [10, 129, 130]). Because of the rapidly increasing 
number of papers on the subject, it is almost impossible to list all the contributions to 
the subject. Therefore, I will just recall the most famous SLs. For example, it has been 
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log(R) Log(R) 
Fig. 3. - R 1 /™ 1 fit. R is in arcsec, /i in mag/arcsec 2 (from [100]). 



found that 

(8) M BH (xC 

where a is the projected central (or within R e ) velocity dispersion of the hosting galaxy; 
after some debate, the currently accepted value of a is very near to 4 (e.g., see [4, 5, 7, 
131]). An important characteristic of this relation (known as the Mbh-c relation) is 
its extremely small scatter, consistent with measurements errors only, so that eq. (8) is 
often considered a "perfect" relation. 

Thus, to a good accuracy, both the Mbh-c and the FJ relations indicate a propor- 
tionality to the fourth power of a , implying the following linear relation between the 
mass of the central SMBH and the total mass in stars M* of the host galaxy 

(9) M BH 1.4 x 1(T 3 M*. 

The latter relation is called the Magorrian relation, after its presentation in [3] , and has 
undergone several successive refinements (e.g., see [9, 132, 133]). Note that versions of 
the above correlation are also provided in terms of the spheroid luminosity instead of the 
spheroid stellar mass [9]. Quite naturally, the possibility of the existence of a FP 
analogous to that of galaxies, but involving Mbh instead of galaxy luminosity, has also 
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Fig. 4. - Panel a: FJ relation for the observed clusters (filled triangles, dashed line) and the dark 
matter halos obtained from cosmological simulations when eq. (84) is used for the mass-to-light 
ratio (empty circles, solid line). L is given in 10 1 L©//i , o in 1000 km/s, R c in Mpc//i, for 
an Hubble constant value of Ho = 100ft km s _1 Mpc -1 , with h — 0.7 (see Sect. 5.3). Panel 
b: with the same symbols, the Kormendy relation for the same data as in Panel a. Each dark 
matter halo is represented by 3 empty circles corresponding to the 3 line-of-sight projections, 
(from [143]). 



been explored, in search of a relation possibly tighter than eqs. (8)-(9). However, so far 
no definite answer has been reached yet (e.g., see [9, 130, 134]-[137]). 

Finally, I wish to mention there is another family of interesting SLs relating the SMBH 
mass to the galaxy Sersic index [138, 139]. For example, [6, 140] found that 

(10) log M BH ^ 2.69 log(m/3) + 7.81, 

with scatter as small as that of the Mbh-c relation. 

2'3. The scaling relations of clusters of galaxies. - The FP is not an exclusive property 
of early type galaxies. For example, clusters of galaxies (within the limitations of a poorer 
statistics) also define their own FP (e.g., see [141]- [143]). In addition, other SLs exist 
for clusters, such as those based on the luminosity and temperature of the Intracluster 
Medium (e.g., see [144]-[147]), but they will not be discussed here. The existence of 
a FP for clusters (for example defined in terms of their total luminosity in stars, of 
the effective radius of the galaxy projected number density, and on some characteristic 
velocity dispersion of the galaxy population) is interesting because the physics behind 
cluster formation is expected to be quite different from the (more complex) physics 
of galaxy formation. Therefore, important clues can be obtained by comparing the 
SLs of clusters and those of galaxies: this will be one of the main subjects of Sect. 5. 
In particular, Lanzoni et al. [143] performed a Principal Component Analysis on the 
Schaeffer sample of Abell clusters, obtaining the following cluster FJ and Kormendy 
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Fig. 5. - The FP of observed clusters (filled triangles) and cosmological dark matter halos when 
eq. (84) is used for the mass-to- light ratio (empty circles) (from [143]). 



relations: 

(11) Loc^- 18 ^ 52 , L oc i?^ 55±0 - 19 , 

where the errors on the exponents also include the observational uncertainties (see 
also [148]). As can be seen from Fig 4, where the relations given in cqs. (11) are plot- 
ted together with the Schaeffer data, the two relations above describe scalings among 
the cluster properties, even if their scatter is quite large (rms dispersion 0.19 in both 
cases). Similarly to what happens for galaxies, a considerable improvement is achieved 
by combining all the three observables together in a FP relation 

(12) Loci?^ ' 15 ^- 31 "- 22 . 

This relation is shown (edge-on) in Fig. 5, and is characterized by an rms dispersion of 
<~ 0.07. Thus, the similarity of the FP for clusters with that for galaxies is remarkable. 
For example, the FP of Es in the B band can be written as L cx i?~ 8 a" 1 - 3 (e.g. [41, 68, 
74, 149]). The situation is similar for the Kormendy relation: in fact, L oc i?J- 7±0 07 has 
been reported for ellipticals in the B band ([79]). The only apparent difference is that of 
the FJ: in fact, Lau 4 for (luminous) Es ([79]), while for clusters the slope is around 2. 

2'4. Cooling flows and AGN underluminosity . - Ellipticals present a well-known "cool- 
ing flow problem" : the time for a significant fraction of the centrally located gas to cool 
(via the observed radiative output) and to collapse is short (typically ~ 10 6 — 10 7 years) 
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compared to the age of these systems (e.g. [150]), while in Es the mass return rate from 
the (passively) evolving stellar population (the main source of Interstellar Medium) is 

(13) M.(t)~1.5x ^ tr 5 13 M yr-\ 

where Lq is the present galaxy blue luminosity and iis is time in units of 15 Gyr ([151]). 
It is then clear that a long-lived cooling flow would accumulate a mass in the galaxy cen- 
tral regions substantially exceeding that currently observed for SMBHs or in the resident 
diffuse gas. Young stellar populations observed in the body of ellipticals also cannot 
account for the total mass released and alternative forms of cold mass disposal (such 
as distributed mass drop-out/star formation) arc not viable solutions (e.g. [152]). In 
addition to this mass disposal problem, which must be solved in any scenario of SMBH 
formation and evolution, the X-ray luminosity L x of low-redshift Es is also inconsistent 
with the standard cooling flow model. In fact, low-redshift Es with optical luminos- 
ity Lb>3 x 10 10 L Q show a significant range in the ratio of gas-to-total mass at fixed 
L B , with values ranging from virtually zero up to few % ([153]); most of the gas is 
detected in X-rays, with temperatures close to the virial temperatures of the systems 
(~ 10 6 - 7 K, e.g. [154]). The idea of a SMBH feedback^) on the Interstellar Medium is 
then most natural, as will be introduced and discussed in Sect. Here I just recall 

that there is increasing evidence in the local Universe of hot gas disturbances on various 
galactic scales, most likely resulting from recent nuclear activity (e.g. [155]). For exam- 
ple, Chandra revealed two symmetric arm-like features across the center of NGC4636 
([156, 157]), accompanied by a temperature increase with respect to the surrounding hot 
Interstellar Medium. The existence of these features has been related to shock heating 
of the Interstellar Medium, caused by a recent nuclear outburst. Other pieces of evi- 
dence include the observation of: a hot filament in the nuclear regions of NGC821 and 
NGC3377 ([158]-[160]); a "bar" feature, presumably due to a shock, at the center of 
NGC4649 ([161]); a nuclear outflow in NGC4438 ([162]); an unusual temperature profile 
present in NGC 3411 ([163]). 

I finally mention another observational riddle that it is almost certainly associated 
with some kind of feedback, i.e., that of the apparent "undcrluminosity" of Active 
Galactic Nuclei in the local Universe, which are observed to emit less than expected, 
based on standard stationary accretion models applied to the circumnuclear environ- 
ment (e.g. [164]- [167]). Quite obviously, the cooling flow and undcrluminosity problems 
should be interpreted and solved in terms of the co-evolution of SMBHs and of the host 
spheroids. In fact, any plausible scenario for the initial formation of the host galaxy and 
of the central SMBH, should also be accompanied by the identification of the mechanism 
able to stop the SMBH growth over the entire galaxy life. 

3. The FP and its interpretation 

In this Section, following [168], I address the problem of the physical interpretation of 
the FP (and of the other SLs) followed by early-type galaxies. This because the galaxy 



( 2 ) The term feedback generally describes the result of the physical phenomena associated with 
the interaction of the emerging radiation from the accreting material onto the host system, on 
different spatial and temporal scales. 
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SLs, when considered in the context of galaxy formation, provide important constraints 
also on the possible evolution of their central SMBHs. 

Quite surprisingly, despite the large amount of dedicated work, no definite interpreta- 
tion of the FP in terms of the intrinsic galaxy properties has been found yet. In general, 
studies of the FP (e.g. [169]-[179]) are carried out under the guiding principle that the 
FP reflects the existence of an underlying mass-luminosity relation for such galaxies 
(e.g. [180, 181]), in a scenario where galaxies are homologous systems in dynamical equi- 
librium. However, as we will see this is not necessarily the case. Given the importance 
of this subject, I review here some of the proposed solutions. I start by giving a short 
discussion about the relation between the FP and the Virial Theorem. Then, I present 
various possible causes to the origin of the observed systematic departure from homology 
(structural, dynamical) and/or from a constant stellar mass-to- light ratio. 

The characteristic dynamical time of Es (e.g., within R e ) is Td yn ~ (G < p > e )~ 1,/2 ~ 
10 8 yrs and their collisionless relaxation time is of the same order ([182]), i.e., both are 
short with respect to the age of Es. As a consequence, presumably only very few galaxies 
currently undergoing strong perturbations are caught in a non-stationary phase, while 
most Es obey the Virial Theorem which, for a galaxy of total stellar mass M* embedded 
in a dark matter halo of total mass Mh, can be written as 

(14) ^^ jq^^ gu.^,.^^ 

where ay and are the so-called (three-dimensional) virial velocity dispersion and total 
kinetic energy of the stellar component, 1Z = M^/M* and T* = M^/L-q is the stellar 
mass-to- light ratio in the specific band used to measure Lb and R e . The dimcnsionless 
functions f/** and W*h are the stellar gravitational self-energy and the interaction energy 
between the stars and the dark matter halo, corresponding to 

(15) £/** = - J < x, V$* > /9*rf 3 x, W*h = -J < x, V$ h > p*d 3 x: 

where $* and <&h are the gravitational potentials generated by stars and dark matter, 
respectively, and <, > is the standard scalar product. Thus, the r.h.s. of eq. (14) depends 
only on the luminous and dark matter density profiles (with the dimensionless functions 
of the order of unity, see, e.g. [102, 183, 184]). Obviously, the virial velocity dispersion 
is a global measure of the total kinetic energy, and does not depend on how such energy 
is distributed between ordered or disordered motions nor on whether the pressure tensor 
is isotropic or not. 

In turn, a v > is related to through a dimcnsionless function that depends on the 
galaxy structure, its specific internal dynamics and on projection effects: 

(16) <7 V = Ck [structure, anisotropy, projection] x a^. 

It is important to note that - even in the case when the pressure tensor of the stellar sys- 
tem is isotropic - Ck is very sensitive to galaxy-to-galaxy structural differences, because 
it relates a weakly structure dependent quantity (cry) to a local property (the central 
projected velocity dispersion a Q ). In fact it can be easily proved that using larger and 
larger apertures to define the relevant projected velocity dispersion <r a (see eq. [27]), in 
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a spherical system without dark matter Ck — ► 3, independently of the galaxy internal 
orbital structure ([185]). If we define 



(17) 

from eqs. (14)-(16) we obtain 



C K 



|t/*J+ft|W* h | 



(18) ^^ = Ky*l 
From eq. (18) and eqs. (2), we get 

(1 9 ) fc i = -7= lo g — ^ — ; fc 3 = -7= log ■ 



while using eq. (1), we find 

Y 2-103 + a 5<3-o 

(20) — ^ oci? e ° Lb~- 

Ay 

iVoie that the Virial Theorem does not imply any FP. In fact, for fixed Lb different 
galaxies, all satisfying the Virial Theorem, can in principle have very different Ky and 
T», and so be scattered everywhere in the fc-space. Equation (20) reveals instead that 
the FP cannot result from pure homology. In practice, in real Es, no matter how complex 
their structure is, the quantity T^/Ky is a well-defined (i.e., with little scatter) function 
of any two of the three observables (L,R c ,cr ). This systematic trend is known as the 
FP tilt (see also Sect. 2.1.1). Note the curious fact that while in the B-band all the 
tilt depends on luminosity (e.g. [72, 100]), in the K-band it is almost due only to i? e , 
with Y*/Ky cx _L^ 02 i?°' 28 . It is evident that fine tuning is required to produce the tilt, 
and yet preserve the tightness of the FP. The central problem posed by the existence of 
the FP is to determine what is the specific physical ingredient responsible for the tilt. 
Of course, the identification of such ingredient would be of extreme importance for our 
understanding of galaxy formation. 

A first possibility to interpret the trend in eq. (20) is to assume homology (i.e., Ky 
identical for all galaxies), a variable stellar mass-to-light ratio (e.g. [181, 186]), and the 
empirically suggested identity 2 — 10(3 + a = 0. Under these assumptions the relation 
would become 

(21) T„cxL B ^, 5 = ^— ^ ~ 0.30 ± 0.064. 

2a 

This line of investigation is addressed in Sect. 3.1. 

An alternative extreme explanation (Sect. 3.2) can be proposed by assuming a con- 
stant stellar mass-to-light ratio and the existence of weak homology ([90, 94, 97], 
[186]-[189]). In this case it is the quantity Ky that is required to be a well-defined 
function of R c and L. If we take 2 — 10/3 + a = 0, the required dependence is 

(22) Ky oc L B - S , 



with S and its scatter the same as above. Equations (21)-(22) represents what Renzini 
& Ciotti ([186]) called "orthogonal explorations". 
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3'1. A stellar origin: changing the Initial Mass Function. - When considering the 
possibility of a stellar origin for the FP tilt, a first candidate is certainly a systematic 
increase of stellar metallicity with galaxy mass. In fact, brighter galaxies are more 
metal rich and then redder than fainter galaxies, so that the stellar mass-to-light ratio 
in the blue band increases systematically with galaxy mass, as requested by eq. (21). 
However, this effect accounts for only a minor fraction of the observed tilt, because the 
FP tilt is still present when using observations at longer wavelengths, such as the K band 
(which is a very good proxy for the almost metallicity-independent bolometric luminosity, 
e.g. [41, 190]). 

A second possible explanation of the FP tilt based on the properties of stellar pop- 
ulations is a systematic change of the so-called Initial Mass Function, i.e. the number 
distribution of stars in a galaxy as a function of their mass; observationally the Initial 
Mass Function of stars in the Solar neighborod is fairly well described by a power law 
of stellar mass (with negative exponent). As is well known from stellar evolution, stellar 
luminosity increases more than linearly with stellar mass, so that the mass-to-light ratio 
of low mass stars is much higher than the mass-to-light ratio of massive stars. Therefore, 
a formal solution for the FP tilt could be a systematic change of slope in the Initial Mass 
Function (with more negative slopes in high-mass galaxies), or a decrease of the minimum 
stellar mass with galaxy mass (while maintainig the Initial Mass Function slope similar in 
all galaxies). These formal solutions were studied in detail in [186]. The main conclusion 
was that a major change of Initial Mass Function slope in the lower main sequence with 
galaxy mass is necessary to account for the FP tilt. At the same time, in order to preserve 
the small and constant thickness of the FP, the galaxy to galaxy dispersion in the Initial 
Mass Function properties should be extremely small, otherwise the edge-on view of the 
FP would rather look like a wedge (with the wide side at high luminosity), rather than 
a strip of constant thickness. Such very small galaxy to galaxy dispersion, coupled to a 
large systematic variation of the slope of the Initial Mass Function, is a rather demanding 
constraint, showing that fine tuning is required to produce the observed tilt of the FP, 
while preserving its constant thickness: the Initial Mass Function should be virtually 
universal for given galaxy mass, and yet exhibit a significant variation with galaxy mass. 

We finally note that, besides of the Initial Mass Function, T* is also a function of age, 
because so is the galaxy luminosity. Thus, the small thickness of the FP implies a small 
dispersion in the age of the bulk stellar content of cluster Es. For example, if galaxies 
are older than ~ 10 Gyr, one gets a tight constraint on the age dispersion, with an upper 
limit of only 1.5 Gyr. Such small age dispersion is in agreement with the similarly tight 
limit that is independently set by the tightness of the color-cr relation for Virgo and 
Coma ellipticals ([42]). 

3 '2. A structural/dynamical origin. - In this case, by assuming =const., one would 
like to explore under which conditions structural/dynamical effects may cause the FP 
tilt via a systematic decrease of Ky with luminosity. Most of these explorations refer to 
spherical, non rotating, two-component galaxy models, where the light profiles resemble 
the i? 1 / 4 law when projected (e.g. [103, 186, 187, 191]-[193]), while very few investi- 
gations consider non-spherical models (e.g. [181, 194, 195]). Bertin in [2] presented a 
lucid discussion of the four main different ways of working on the construction of stellar 
dynamical models. Here I will focus on the so-called p — to — f or "density priority" 
approach. In the spherical cases, usually carried out under the assumption of a two- 
integral phase-space distribution function / = f(E,J 2 ) (with E and J being the star 
specific energy and angular momentum; in this case the tangential components of the 
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velocity dispersion tensor are identical, the only possible difference being between of and 
a e = °0 = °f/2' c -g- [196]), the model spatial and projected velocity dispersion profiles 
are obtained by solving the associated Jeans equation for assigned density and pressure 
anisotropy profiles (see, e.g. [197]): 

dp*(r)<j 2 (r) 2[3(r)p*(r)<ji(r) GM(r) 
[26) Jr + r ~ P * [ )l 

with the boundary condition p*(r)<7 2 (r) — > for r — > oo, where M(r) is the total (dark 
and visibile) mass within r. Usually the orbital anisotropy is introduced in terms of the 
Osipkov-Merritt ([198]- [200]) parametrization 

( 24 ) 0{r) = 1 - 4?4 = -T^' 

o?(r) r 2 +r 2 J 

so that the velocity dispersion tensor is nearly isotropic inside r a and radially anisotropic 
outside. The radial component of the velocity dispersion is then given by 

(25) P*(rK 2 M = ^2 J" P*(r)M(r) (l + if) dr, 
The projected velocity dispersion is given by 

(26) ^ = ^ /, 



Vr 2 - R 2 



where is the surface stellar mass density. The "observed" a a does not corre- 

spond to <7p(0), but rather to the average over the aperture used for the spectrographic 
observations: 

(27) v 2 M= j^^J^Z*(R)o$(R)RdR, 

where Mp(R a ) is the projected stellar mass inside R a . Usually, one refers to i? a = Re/8 
or fi e /10. 

Examples of density profiles that have been extensively assigned are those obtained 
by deprojecting the i? 1 /™ profiles (as defined in eq. [6]), and the so-called 7- models 
([184, 201]) 

(28) *( r ) = ^_J^ (0< 7 <3), 

which for 7 = 1 (R c ~ 1.82 r*) and 7 = 2 (i? c ~ 0.75 r*) reduce to the well known 
Hcrnquist ([202]) and Jaffe ([203]) profiles, respectively. In order to study the effects of 
dark matter halos, dark matter can be introduced to contribute to the total integrated 
mass M(r) and be represented again by 7 models (with different total mass and scale- 
length [191, 193]), or by the Plummer ([204]) density profile 



CO-EVOLUTION OF ELLIPTICAL GALAXIES AND THEIR. CENTRAL BLACK HOLES 



15 



and its variants, such as the quasi-isothermal density profile p oc (r 2 + r 2 ) 1 and the 

King ([205]) p oc (r£ + r 2 ) density profile. 

In general, two-component 7 models are associated with a density profile in the central 
regions similar to that expected from cosmological simulations, i.e, with a central cusp 
(e.g. [206]- [209]). In addition, recently two-component models in which the total mass 
profile decreases as r~ 2 have been found consistent with combined observations of stellar 
dynamics and gravitational lensing (e.g. [210, 211]). 

It should be remarked that not all the models adopted in the FP investigations have 
realistic velocity dispersion profiles in the central regions. Many of them (for example, 
those based on Hcrnquist or R 1 /" 1 density profiles) from the solution of eq. (23) are 
characterized by a sizable central depression in their erp (for a detailed discussion of this 
point, see [94, 100, 102, 103, 187, 212]), while the observed profiles typically decrease 
monotonically with radius (e.g. [213]-[217]). Thus, some care is needed in the model 
selection. 

Among models with realistic velocity dispersion profiles are the isotropic Jaffe models 
and the /oo models ([61, 218], and references therein). In particular, the /oo models have 
been constructed by following the physical scenario that Es may have formed through 
collisionless collapse ([60]). In the spherical limit, their anisotropic distribution function 
is given by 

, _ ( A(-E) 3 / 2 cxp(-aE-cJ 2 /2) if E < 0, 
(M) /oo "\0 if£>0, 

where E = v 2 /2 + $(r) and J are the star energy and the magnitude of the star an- 
gular momentum, per unit mass, respectively; here A, a, and c are positive constants. 
The models are a one-parameter family, with structure dependent on the concentration 
parameter 5» = — a$(0). 

3'2.1. Dark matter content and distribution. In this line of investigation, in the con- 
struction of the galaxy models one usually assumes isotropic velocity dispersion, and 
ascribes all the FP tilt to a systematic variation with galaxy luminosity of the relative 
amount of dark matter with respect to the stellar component of the galaxy (parametrized 
by 1Z = Mh/M*) or to a systematic variation with galaxy luminosity of the relative con- 
centration of the dark matter (parametrized by — r^fr*, where is a characteristic 
radius of the dark matter distribution). In the case of models with fixed 0, the larger 0, 
the larger the variations of 1Z that are required to produce the tilt (see the top panels in 
Fig. 6 for two examples). In any case, in this interpretation of the FP tilt massive galax- 
ies should be more dark matter dominated than smaller systems. In the complementary 
exploration, one can fix the value of the total mass ratio 1Z = M^/M* and explore the 
required variations of = r^/r^: in general, it turns out that should decrease as a 
function of galaxy mass (bottom panels in Fig. 6). 

Overall, as already stressed, the narrow and nearly constant thickness of the galaxies 
distribution about the FP corresponds to a very small dispersion of T^/Ky, and if T* and 
Ky are not finely anticorrelated, this implies indeed a very small dispersion, separately 
for both quantities, at any location on the FP. In the present context (T* =const), this 
translates into strong constraints on the range that 1Z and can span at any location on 
the FP. For example, the dotted lines in Fig. 6 represent the band within which galaxy 
to galaxy variations of the corresponding parameter are allowed, and yet are consistent 
with the restrictions imposed by the tightness of the FP. It is evident from these figures 
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Fig. 6. - The trend along the FP of the dark matter content (upper panels) at constant (3 and 
that of the dark matter concentration (lower panels) at constant 1Z, required to produce the tilt, 
in two-component isotropic Hernquist models (HH) and two-component isotropic Jaffe models 
(JJ) isotropic models. The band within dotted lines marks the boundaries within which 1Z 
and (5 can vary at each location on the fci axis in accordance with the observed FP tightness 
(from [94]). 



that, whatever the structural parameter that is responsible for the tilt of the FP, and 
whatever the assumed mass distribution, dramatic fine tuning is required to produce the 
tilt, and yet preserve the tightness of the FP. 

3'2.2. Structural non-homology. A change of internal structure as systematic way to 
change Ky has been proposed and studied by several authors (e.g. [86, 90, 94, 96, 97, 100, 
103, 186, 187, 188, 219, 220]). We may call weak homology the condition by which the 
structure and dynamics (density and pressure tensor distributions) of early-type galaxies 
vary systematically with galaxy luminosity. 

An indication that this may be at the origin of the FP tilt is indeed provided by 
the systematic trends noted above in terms of i? 1 /" 1 fits (see Sect. 2.1.3). Therefore, 
this option can be investigated using isotropic i? 1 /" 1 models without dark matter ([102]), 
thus ascribing the origin of the tilt to a systematic variation of m. For example, in the 
preliminary analysis of [94] it has been shown that in order to produce the tilt m has to 
increase from 4 up to <~ 10. These values are well within the range of values spanned by 
observations. 
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Fig. 7. - The values of m for the six f? 1//m isotropic models required to produce the FP tilt, and 
the band within which m can vary at each location on the fei axis in accordance with the observed 
FP tightness. In the upper panels, the faintest model (lowest value of k\) is characterized by 
m — 4; in the lower panel, m = 2 at the faint end of the FP (from [94]). 



A much more in-depth analysis of weak homology has been carried out in [100], 
starting with a very accurate analysis of 4 significant cases of deviations from the i? 1 / 4 
law, i.e. the surface brightness profiles of NGC 1379, NGC 4374, NGC 4458, and NGC 
4552. The analysis (based on the photometric data taken from [112, 113]) confirms 
that the i? 1/<m can provide superior fits (the best-fit value of m can be lower than 2.5 
or higher than 10), better than those that can be obtained by a pure i? 1 / 4 law, by 
an .R^+exponential model, and by other dynamically justified self-consistent models 
(such as the models ([218]). Note that according to eq. (22) a factor of 20 in Lb 
would require a change in K\ by a factor of w 2.45. Then, it is of interest to compute 
the function Ky(m) for one-component, spherical, non-rotating, isotropic i? 1 /" 1 models 
and a simulated spectroscopic aperture of R a /8 [103, 97]. An accurate and convenient 
analytical representation in the range 1 < m < 10 (with typical errors on the order of 
one percent) is given by 

73 32 

< 31 > *v(m)* 10465 +(m _ . 94)2 + 0-954- 

A similar argument also holds for the models, for which the simple analytical inter- 
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polation formula 

142.3 - 41.51* + 2. 66* 2 



(32) ifv(t) 



30.61 - 10. 7* + * 2 



is accurate to around one percent, in the range 2 < \& < 10. Incidentally, this shows 
again how the determination of the coefficient Ky is sensitive to the choice of models 
(e.g., see Fig. 11 in [100]). Note that this leads to significantly different answers with 
respect to the application of R}/ m models. For example, if we take NGC 4552, the R 1 /™ 1 
modeling would give m ~ 11.14 and Ky ~ 1.6, while the /oo modeling would set * ~ 6.2 
and Ky ~ 2.5; for the galaxy NGC 4458 (closest to the standard i? 1 / 4 law in our sample), 
we find to ~ 3.2 and Ky ~ 5.7 based on i? 1//m modeling or \P ~ 9, Ky ~ 3 based on 
modeling. In any case the variations found are consistent with the FP tilt. 

A more sophisticated analysis of the problem, based on Monte-Carlo simulations, was 
also presented in [100]. In practice, the Authors mapped the model space (to, T„, L, R e ) 
into the observed space (L,R c ,cr ) by using the virial coefficients in cqs. (31)-(32), and 
selecting as acceptable candidates for real galaxies, only those points that turn out to 
be compatible, in the observed space, with the FP correlation and its observed scatter. 
For example, Fig. 8 shows, under the assumption of constant T„, the distribution in 
the luminosity-Sersic index (m) plane of galaxies compatible with the FP: in particular, 
models represented by heavy dots are consistent with a FP of virtually no scatter (what 
we may call the "backbone" manifold). Of course, by allowing a simultaneous variation 
of T* and m, a much larger number of models becomes compatible with the FP. This case 
is represented in Fig. 9, where all the points are compatible with the FP with intrinsic 
scatter artificially reduced. Thus, these results show that weak homology - both from 
the observational and from the theoretical point of view - could be a relevant physical 
ingredient at the origin of the FP. 

3'2.3. The role of anisotropy. Among the various galaxy properties in principle able to 
destroy the FP thinness (as a consequence of a substantial variation at fixed galaxy lumi- 
nosity), one possibly "effective" might be orbital anisotropy (e.g. [221]). In fact, galaxy 
models are often believed to be able to sustain a large spread of orbital anisotropies. It 
is also well known that radial orbital anisotropy can be associated with very high central 
velocity dispersion values, and correspondingly low values of Ky, thus violating the FP 
thinness in case of significant scatter among galaxies. Natural questions to be addressed 
are then 1) if varations in orbital anisotropy can be responsible for the whole FP tilt, 
and 2) whether empirically, or as a result of some physical reasons, the range of orbital 
anisotropies present in real galaxies might be limited. In practice, the amount of radial 
anisotropy in the velocity dispersion tensor should increase with galaxy luminosity (i.e., 
r a should decrease in the Osipkov-Merritt parameterization), but with values fine-tuned 
with galaxy luminosity. In other words, in this scenario the FP tilt would be produced 
by a dynamical non-homology due to anisotropy. 

The results of a preliminary investigation ([94]) of this problem, based on the be- 
havior of Hernquist and Jaffc density profiles, are shown in Fig. 10. The curves are 
truncated because of the limits imposed by dynamical consistency ([187, 191, 193]): 
above a certain luminosity, in models constrained to the FP, the phase-space distri- 
bution function would run into negative values. A more refined exploration, based on 
one-component i? 1 /™ models, constructed with different amounts of pressure anisotropy, 
following the Osipkov-Merritt prescription can be found in ([103]). For these models the 
self-consistently generated phase-space distribution function has been obtained, and the 
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Fig. 8. - Projection on the plane (L,m) of the manifold in model space identified by T, ~ 3.5. 
The scatter of points reflects the adopted scatter around the FP in the observed space. Solid 
symbols represent the "backbone" manifold. Open squares represent the four observed galaxies 
described in this Section (from [100]). 



minimum value of the anisotropy radius for the model dynamical consistency has been 
derived as a function of to. As for Hernquist and Jaffe models, also for i? 1 /" 1 models 
constrained to the FP it is found that above a certain luminosity (dependent on to), the 
phase-space distribution function runs into negative values (as shown in Fig. 11 for the 
case m = 4). Therefore, anisotropy alone cannot be at the origin of the tilt, because 
the extreme values of r a that would be required correspond to dynamically inconsistent 
models. 

The relation between radial anisotropy and FP thickness was also studied in [103] in 
a semi-quantitative way by considering the radial-orbit instability indicator £ = 2K r /K t 
([222]), where K t = 2n J pa 2 r 2 dr and K T = 2ir J pa 2 r 2 dr are the tangential and the 
radial kinetic energies (see also [223]). In particular, it is empirically known that when 
£>1.5 -j- 2 a radially anisotropic spherical system is likely to be unstable( 3 ), and in [103] 
was argued that i? 1 /™ galaxy models sufficiently anisotropic to be outside the FP observed 



( 3 ) Unfortunately, such indicator is not fully reliable, because indications exist that it can 
depend significantly on the particular density profile of the model under investigation (e.g. [224]- 
[229]. See also the recent study [230]). 
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Fig. 9. - Projection on the (L, m) plane of the entire backbone manifold identified by the FP 
with scatter artificially reduced to O.Of. Open circles correspond to models with 3.4 < T* < 3.6, 
open squares as in Fig. 9 (from [100]). 



thickness (when their parent isotropic model was assumed to lie on the FP) should be 
unstable. The relevance of this result in connection with the FP thickness is simple: the 
effect on the projected velocity dispersion due to the maximum orbital anisotropy allowed 
by the stability requirement is well within the FP thickness, and so no fine-tuning for 
anisotropy is required. 

Of course, those above were qualitative expectations based on necessarily simple an- 
alytical models. The entire question was then addressed quantitatively by using high 
resolution A^-body simulations ([231]) of radially anisotropic one-component and two- 
component 7 models, and exploring the impact of radial orbital anisotropy and instabil- 
ity on the FP properties. The numerical results confirmed the previous studies, and the 
situation is summarized in Fig. 12. The globally isotropic parent models are obtained 
fixing their structural properties and assuming isotropic pressure, and are placed on the 
FP by assigning the pair (T*,Lb): with the coordinates adopted in the figure, the par- 
ent models are placed on the origin. From each of these parent models lying on the FP 
one then generate a family of Osipkov-Merritt radially anisotropic initial conditions by 
decreasing s a = r a /R e , while keeping all the other model parameters fixed. The dotted 
line is the locus of points corresponding to radially anisotropic initial conditions, and for 
sufficiently small values of s a the members of each family are found outside the observed 
thickness of the FP. The quantity in the vertical axis corresponds to the distance from 
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Fig. 10. - The trend of the anisotropy radius in units of R c along the FP required to produce 
its tilt, in Hernquist and Jaffe models. The curves are truncated at the value of r a /R e below 
which the models become dynamically inconsistent. The band within dotted lines marks the 
boundaries within which r a can vary at each location on the k\ axis according to the observed 
FP tightness (from [94]). 



the FP of each virialized end-product of the numerical simulations: if k% = kf° then the 
model has "fallen back" on the FP, so points inside this strip represent models consistent 
with the observed thickness of the FP. It is apparent how, by increasing the amount 
of radial anisotropy, the initial conditions move along the dotted line, and when they 
reach the critical value of the anisotropy radius they become unstable( 4 ) and rearrange 
their density profile and internal dynamics in a new, stable configuration. It is apparent 
that all the models, with the exception of two, obey the predictions made in [103]: the 
restrictions set by the onset of the radial-orbit instability match almost exactly the FP 
thickness. 

Moving to the strictly related question whether the FP tilt can be due to a systematic 
variation of Ky induced by an appropriate underlying correlation sa-Le at fixed galaxy 



( 4 ) In [231] it is found that for one-component 7-models the critical value £ s for the radial-orbit 
instability is in the range 1.6<£ s <1.8. This range for £ s is compatible with the value 1.7, 
reported by [222] and used in [103], and with the results of [228], who estimated £ s ~ 1.58 for 
the family of f x models; [229] found instead a higher threshold value for stability (£ a ~ 2.3). 



22 



L. CIOTTI 




Fig. 11. - The modifications of the (dimensionless) R 1 ^ 4 distribution function as a function of 
the Osipkov-Merritt parameter Q — E + J 2 /(2r a ) (normalized to the system central potential), 
moving from the globally isotropic case (solid line) to the critical anisotropy for consistency 
(dashed line), where s a = r a ./R e . This behaviour is common to the whole family of the B}l m 
models, and seems to be more a general property of the adopted anisotropy profile rather than 
a characteristic of some specific mass model (see also [191, 193]) (from [103]). 



structure, in [231] it is also shown that a systematic increase of radial orbit anisotropy 
with galaxy luminosity is not a viable solution, because the galaxy models become unsta- 
ble at moderately high luminosities, and the virialized systems originated from unstable 
initial conditions fall well outside the FP itself. This is illustrated in Fig. 13, where 
k\ — kf° measures how much a given initial condition is displaced on the FP from its 
parent isotropic model, and fc 3 — fcg measures the position of the corresponding end- 
product. It is clearly impossible to reproduce the FP tilt over the whole observed range 
(e.g., Afci ~ 2 and Ak^ ~ 0.3 as reported in [75]) by using stable models only. More- 
over, Fig. 13 shows also that the end-products of unstable initial conditions fall well 
outside the FP thickness, and the departure is larger for larger distance from the parent 
galaxy. Thus, we can safely conclude that the FP tilt cannot be explained as an effect of 
a systematic increase of radial anisotropy with galaxy luminosity under the assumption 
of structural homology. 

3'3. Summary. - In this Section the main plausible hypotheses proposed to interpret 
the properties of the FP have been briefly illustrated. In general, it is shown that, 
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Fig. 12. - Final vs. initial k-$ for one- (solid dots) and two-component (empty dots) galaxy 
models obtained from TV-body simulations, where kf° is the coordinate of their isotropic parent 
galaxy. The horizontal dashed lines mark the FP observed thickness (j(kz), while the dotted line 
&3 = k\ is the locus of the initial conditions. The end-products are asymmetric, and so their 
representative points span a range of values as a function on the line-of-sight orientation (vertical 
lines); remarkably, the length of these segments is smaller than the FP thickness (from [231]). 



independently of the favoured explanation, a remarkable fine-tuning of galaxy structure, 
dynamics, and properties of stellar population (in particular, of the Initial Mass Function) 
with galaxy luminosity must exist in order to produce the FP tilt and yet to preserve 
its thinness. Overall, the most promising explanation of the FP tilt seems to be 1) a 
synchronized star formation of the bulk of stars; 2) a remarkable regularity in the amount 
and distribution of dark matter; 3) a systematic dependence of the galaxy stellar density 
profile on galaxy luminosity (weak homology). This conclusion has been reached also 
in [232]. Anisotropy seems to be ruled out as a significant contributor to the FP tilt, as 
well as to its thickness. This latter conclusion is reached both on theoretical grounds (i.e., 
stability arguments) and from powerful empirical evidence, i.e., the very existence of the 
Mbh-Co relation. In fact, no significant scatter in the anisotropy can be present in real 
Es, otherwise the thin Mbk-cto relation would be destroyed (by the way, the thinness of 
this latter relation also shows that projection effects, due to non sphericity of galaxies, 
cannot affect significantly the observed a Q , see also [194]). 



24 



L. CIOTTI 



n — i — i — i — I — i — i 1 — i — I — i — i — i — i — I — i — i — i — i 1 — r 



0.05 



■» * ■ J ",fo 



-0.05 

■ i — i CO 

I 

CO 

^-0.1 



-0.15 



-0.2 



_i i i i i i i i i i i i i i i i i i i i i_ 



0.05 



0.1 



0.15 



0.2 



Fig. 13. - End-products of iV-body simulations with anisotropic initial conditions placed on the 
FP at kg — kf°, where kf is the coordinate of their isotropic parent galaxy. The FP observed 
thickness a $3) is marked by the dashed lines (from [231]). 



4. Galaxies and SMBHs: maintenance of Scaling Laws 

In this Section two important phenomena, in principle able to destroy the galaxy 
and SMBHs scaling laws, are discussed. The first, galaxy merging, is of external origin, 
while the second, namely the disposal of the huge (compared to the observed mass of the 
SMBHs) mass return from the passively evolving stellar population in Es, is of internal 
origin. 



4T. Merging. - One of the basic problems posed by dry merging (i.e. non-dissipative 
merging of red, old galaxies with no significant amount of gas) as the main channel to 
form Es was explicitly pointed out and discussed heuristically in [31]. In practice, it is the 
impossibility to increase the velocity dispersion of the end-products (as requested by the 
FJ law) while, at the same time, the associated enormous growth of R e leads to violate 
the Kormendy relation. The two facts above are direct consequences of basic physics. In 
fact, from the virial theorem and conservation of the total energy, in the merging of two 
galaxies on a parabolic orbit (with masses Mi and M 2 , virial velocity dispersions ov,i 
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and cry, 2, virial radii rv,i and ?"v,2), the virial velocity dispersion and the virial radius( 5 ) 
of the resulting galaxy, in the case of no mass loss and negligible kinetic and interaction 
energies of the galaxy pair when compared to their internal energies, are given by 



It follows that 0V,i+2< max(oy .1, oy^) an d rv,i+2> mm (ry 1, ry 2 ), i.e., the virial veloc- 
ity dispersion cannot increase and the virial radius cannot decrease in a merging process 
of the kind described above. For example, in a merging hierarchy of identical seed galax- 
ies characterized by cry ; o, ry t o and Mo, we expect cry = cy,o and ry = (M/Mo)rv,o ; 
independently of the merging sequence: if a <~ oy and R c ~ ry, then it results that the 
FJ and Kormendy cannot be consistent with dry merging. 

Several high-resolution N-body simulations of dry merging are nowadays available 
(e.g. [219], [233]-[242]) that can be compared with the expected relations (33). Here I 
describe in some detail the results obtained in [236] which, by using numerical simulations 
based on one and two-component isotropic Hernquist galaxy models ([191, 193, 202], see 
also Sect. 3.2), checked whether the end-products of merging of galaxies, initially lying 
on the FP, lie on the FP. 

In particular, the study explored two extreme situations, namely the case of major 
merging, in which equal mass galaxies are involved at each step of the hierarchy, and 
the case of accretion, in which a massive galaxy increases its mass by incorporating 
smaller galaxies. In practice, the first generation of the equal mass merging hierarchy is 
obtained by merging a pair of identical, spherically symmetric and isotropic Hernquist 
models (the "zeroth order" seed galaxies), while the successive generations (in general 
5, for a total mass increase of 32) are obtained by merging pairs of identical systems 
obtained by duplicating the end-product of the previous step. Cases with dark matter 
halos (R = 5 and (3 — 2 in the seed galaxies; for the notation, see Sect. 3.2.1) and non- 
zero orbital angular momentum have also been considered. In the accretion hierarchy, 
instead, the seed test galaxy grows by accretion of smaller systems: the first merging 
event is identical to that in the equal mass merging case, but in the successive steps the 
end-product merges again with a seed galaxy, and so on, until the same final mass as in 
the equal-mass merging case is reached. 

The first important result - which is almost independent of the different cases explored 
- is shown in Fig. 14: while cry is almost constant, the linear growth of ry (and of the 
volume hal-mass radius tm) with M (dotted line) is apparent. The small increase of 
ay with respect to the expectation cry = const, is explained by the modest mass loss 
during merging. It is found that the structure of the end-products is quite sensitive to 
the different growth assumptions. In fact, from Fig. 15 demonstrates that the Sersic 
index m changes with the mass of the merging end-products: while in the equal-mass 
mergers m increases with mass and spans the range 2 < m < 11, as in real galaxies (see 
Sect. 2.1.3), in the accretion case m decreases with the galaxy mass at mass ratios >4 
for the head-on accretions. Therefore, the explored head-on accretion scenario fails at 
reproducing the relation between the surface brightness profile shapes and luminosity of 



( 5 ) By definition, in a one-component galaxy ct v = 2T/M and ry = —GM 2 /U, where T and 
U are the total kinetic and the gravitational energy of the galaxy, respectively. See beginning 
of Sect. 3 for the two-component case. 
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Fig. 14. - Tbp panel: virial velocity dispersion of the stellar component at stage i of the merging 
hierarchy vs. the total stellar mass of the merger. Bottom panel: angle-averaged half-mass 
radius rM vs. total stellar mass. Equal mass mergers are shown as solid triangles and squares 
(one-component galaxies), and stars (two-component galaxies); empty triangles and circles 
represent the accretion hierarchies. Triangles correspond to simulations with non zero orbital 
angular momentum. The dotted line indicates the relation rn oc M; M* ; o, vm,o and <tv,o are 
the stellar mass, the virial velocity dispersion and the half-mass radius of the seed galaxy. In 
the two-component cases, cry and tm refer to the stellar component only. Note the different 
range spanned in the ordinate axes in the two panels (from [236]). 



real galaxies. 

The impact of dry-merging on the FJ and Kormendy relations is summarized in Fig. 16: 
by analogy with the results shown in Fig. 14, the FJ and the Kormendy relations are 
again shown to be violated at high masses (luminosities). Figure 17 shows the position 
of the end-products in the edge-on (k\, fc 3 ) and face-on (k\, fc 2 ) projections of the FP. 
The progenitor of the merging hierarchy (the black dot without bar) is placed exactly on 
the edge-on FP at ki<3. Consistent with the adopted dissipationless scenario, the value 
of T* is kept constant during the whole merging hierarchy. Due to the loss of spherical 
symmetry of the end-products of the merging simulations, their coordinates depend on 
the line-of-sight direction; however, being the luminosity (mass) of each end-product 
fixed, variations of k\, k 2 , k 3 due to projection effects are not independent. In fact, the 
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Fig. 15. - Sersic best-fit parameter m vs. total stellar mass of the end-products at stage i of 
the merging hierarchy. Same symbols as in Fig. 14 (from [236]). 



ki and fc 3 coordinate of a galaxy of given luminosity are linearly dependent as 




corresponding to the segments in the top panel. The main conclusion of this analysis is 
that one-component (solid squares and triangles) and two-component (stars) equal mass 
mergers behave almost in the same way: models climb over the edge-on FP and remain 
well inside the populated zone in the face-on FP, moving along roughly parallel to the 
line defining the zone of avoidance (see also [75]). Thus, equal mass dry merging seems 
to be surprisingly consistent with the existence of the FP, especially considering its small 
thickness when seen edge-on. The most striking difference of the end-products of head- 
on accretion simulations (empty circles) with respect to the equal mass merging hierarchy 
and also to accretion simulations with angular momentum (empty triangles) shows up in 
the edge-on FP: after few accretion events, the end-products are characterized by a k 3 
decreasing for increasing k\ , and the models corresponding to an effective mass increase 
of a factor <~ 12 are found at a distance Sk 3 larger than the FP scatter. This result is 
not surprising. In fact, being the coordinate k 3 a measure of non-homology (for galaxies 
with constant T„), the decrease of the Sersic parameter m with increasing mass reflects 
directly in the unrealistic trend in the (k±, k 3 ) plane. In contrast, in the face-on (ki, k^) 
plane, the end-products of accretion (both head-on and with angular momentum) evolve 
along the same direction followed by equal mass mergers, but with a smaller excursion 
in k\ and k^. Thus, while in the equal-mass mergers the deviations from the FJ and 
Kormendy relations curiously compensate to reproduce the edge-on FP, in the accretion 
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log M,/M, 

Fig. 16. - Top panel: stellar central velocity dispersion (normalized to that of the first progenitor) 
vs. total stellar mass. Points correspond to angle-averaged values, bars indicate the range 
spanned by projection effects. The solid and the two dashed lines represent the FJ relation 
Lb oc at and its scatter, while the dotted line represents Lb oc a't b . Bottom panel: stellar 
effective radius (normalized to that of the first progenitor) vs. total stellar mass. Points and 
bars have the same meaning as in the top panel. The solid line represents the adopted "fiducial" 
Kormendy relation. Note the different range spanned in the ordinate axes in the two panels 
(from [236]). 



case there is not enough compensation, and the FP tilt is not reproduced. 

Finally, we can move to consider the effects of dry merging on the Mbh-ct relation( 6 ). 
It is assumed that each seed galaxy contains a BH of mass Mbh,o, and that each merging 
end-product contains a BH obtained by the merging of the BHs of the progenitors. 
Unfortunately, BH merging is still a poorly understood physical process, in particular 
with respect to the amount of emitted gravitational waves (e.g. [31, 244, 245]), and 
for this reason we consider two extreme situations: the case of classical combination of 
masses (Mbh\i+2 = A^bh,i+^bh,2, with no emission of gravitational waves) and the case 
of maximally efficient radiative merging (Mg H 1+2 = Mg H 1 + Mg H 2 , corresponding to 
entropy conservation in a merging of two non rotating BHs [246, 247]). Figure 18 shows 
the central velocity dispersion of the mergers versus the mass of their central BH in the 
case of classical (panel a) and maximally radiative (panel b) BH merging. As expected 
from the similarity between the FJ and the Mbh-c relations, in the classical case, both 
equal mass mergers and accretion mergers are unable to reproduce the observed relation. 
Thus, as for the FJ, the reason of the failure of dissipationless merging at reproducing 



( 6 ) The simulations of [236] do not consider the central SMBH. However, it has been shown [243] 
that the formation of a binary BH does not modify significantly a a . 
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Fig. 17. - Top panel: the merging end-products in the (ki,ks) plane, where the solid line 
represents the FP relation as given by equation (4) with its observed 1-a dispersion (dashed 
lines). Bars show the amount of projection effects. Bottom panel: the merging end-products 
in the (^1,^2) plane, where the dashed lines define the region populated by real (as given 
in [75]). Each model is represented by a set of points corresponding to several random projections 
(from [236]). 



the Mbh-Co relation is that the end-products are characterized by a too low a a for 
given Mbh, i-e., -Mbh is too high for the resulting a Q . A promising solution to this 
problem could be the emission of some fraction of Mbh as gravitational waves. In fact, 
by assuming maximally efficient radiative BH merging, the points are found remarkably 
closer to the observed relation, even if it appears that the slope of the Mbh-c relation is 
not well reproduced by the end-products of head-on accretion (empty circles). While in 
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Fig. 18. - Panel a: galactic central velocity dispersion vs. BH mass for classical BH merging; ao,o 
and Mbh,o are the central velocity dispersion and BH mass of the first progenitor, respectively. 
The points correspond to the mean value over the solid angle, while the bars indicate the range 
spanned by projection effects. Solid and dashed lines represent the Mbh-o" relation for a = 4 
and q = 5, respectively, while vertical dotted lines show the observed scatter around these 
best-fit relations. Panel b: same data as in panel a, but for maximally radiative BH merging. 
Panel c: stellar mass vs. BH mass for classical BH merging; M*^ is the stellar mass of the first 
progenitor and the solid line represents the Magorrian (Mbh oc Mb u i ge ) relation. Panel d: same 
data as in panel c, but for maximally radiative BH merging (from [236]). 



the classical scenario the Magorrian relation is (obviously) nicely reproduced, in the case 
of substantial emission of gravitational waves the relation between BH mass and bulge 
mass is not reproduced. All the results presented hold under the strong assumption that 
the BHs of the merging galaxies are retained by the end-products, but there are at least 
two basic mechanisms that could be effective in expelling the central BHs. The first is 
related to the general instability of three body systems: if a third galaxy is accreted by 
the end-product of a previous merging before the binary BH at its center merged in a 
single BH, then the escape of the smallest BH is likely. It is clear that if this process 
happens more than a few times, then the Magorrian relation will not be preserved at 
the end (for a detailed discussion of this problem see, e.g. [243], [248]- [250]). A second 
physical mechanism, which could be even more effective in expelling the resulting BH 
from the center of a galaxy merger, is related to the possibility of anisotropic emission 
of gravitational waves. This process, commonly known as the "kick velocity", is directly 
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related to the fraction of BH mass emitted anisotropically during BH coalescence. In 
fact, gravitational waves in fact travel at the speed of light, and so even the anisotropic 
emission of a few thousandths of the mass of the BH binary will produce a recoil (due to 
linear momentum conservation) of the resulting BH with a characteristic velocity of the 
order of, or higher than, the escape velocity typical of massive galaxies. In conclusion, 
it is not obvious that in each galaxy merging the resulting BH will remain at the center 
of the galaxy. 

I illustrate now the results of a Monte-Carlo investigation of merging effects, in which 
also gas dissipation is heuristically taken into account [251]. This is done with the aid of a 
simple yet robust approach designed to model the influence of merging on galaxy structure 
and the consequent effects on the scaling laws followed by Es. In fact, simple physical 
arguments show that gas dissipation should be able to mitigate the problems posed 
by dry merging to the explanation of the observed scaling laws (e.g. [31], [252]- [255]). 
Unfortunately, numerical simulations with gas dissipation are considerably more complex 
than pure N-body simulations (e.g. [253, 254], [256]-[262]), so that simple arguments as 
those following can be of help. 

For simplicity, each elliptical is modeled as a non rotating, isotropic and spherically 
symmetric one-component virialized system (i.e., dark matter is distributed proportion- 
ally to the visible matter), characterized by stellar mass M*, gas mass M gas = aM„ and 
SMBH mass Mbh = 1-tM*; from the observations ([3]) we set p ~ 10~ 3 in spheroids of 
the nearby universe (z = 0). The total energy of a galaxy is then given by 

(35) E = K* + U g + W, 
where 

(36) tf, = lf P* °l d 3 x, U g = ^J Pg Td 3 x , W=\j \p,+ p s )(<£* + %) d 3 x 

are the stellar kinetic energy, the gas internal energy, and the total gravitational energy 
of stars and gas. Here <j*, fee, T, and (m) are the stellar 1-dimensional velocity disper- 
sion, the Boltzmann constant, the gas temperature, and the gas mean molecular mass, 
respectively. 

Under the simplifying assumption that the gas is spatially distributed as the stars 
(i.e., p g = ap*) and that the gas is in equilibrium in the total gravitational field, we have 
$ g = a $*; from the Jeans and the hydrostatic equations we find T = (to)o^/A;b, so that 

(37) W= (1 + afW*, U g = aK*, E=-(l + a)K* = ^ + 0t ^ W., 

where W* is the self- gravitational energy of the stellar component. Finally, the charac- 
teristic one-dimensional stellar velocity dispersion cry and the characteristic radius ry, 
defined as 

(38) tf, = jM„4, |W,| = ^— i-, 

2 r v 

are related (in the case of isotropic R}l m models with 2<m<12) to the observables R e 
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and (j (with spectroscopic aperture Re/ 8): 

r v _ 250.26 + 7.15m a Q ^ 24.31 + 1.91m + m 2 
^ ^ ^ ~ 77.73 + m 2 ' ^ ~ 44.23 + 0.025m + 0.99m 2 ' 

For simplicity, in the following we consider only the case of the parabolic merging of 
two galaxies, so that the total energy of the system is the sum of the internal potential 
and kinetic energies of the two progenitor galaxies; we also assume that no mass is lost 
in the process. During merging, as a consequence of gas dissipation, a fraction r\ of the 
available gas mass is converted into stars. The stellar mass balance equation is 

(40) M* = M*i + M* 2 + v(M g i + M g2 ). 

Furthermore, a new SMBH forms by the coalescence of the two central BHs and a fraction 
frj of the available gas is accreted on it, leading to a BH of final mass 

(41) M BH = (M p BH1 + M BH2 )^ + f v (M gl + M g2 ). 

As in [31, 236], p — 1 corresponds to the classical merging case (no gravitational radi- 
ation), while p = 2 to the maximally radiative case for non-rotating BHs. In equation 

(41) it is implicitly assumed that first Mbh,i and Mbh,2 merge, and then the gas is 
accreted on the new BH; the other extreme case would be that of gas accretion followed 
by merging (e.g. [263]). As a consequence of star formation and BH accretion, the gas 
mass balance equation is 

(42) M gas = (l- V - f V )(M gl +M g2 ), 

which implies that < 77 < 1/(1 + ./*)■ Thus, the gas-to-star mass ratio after the merger 
and the new Magorrian coefficient are given by 

(43) a = M gas _ (1 - n - /r?)(aiM„i + a 2 M* 2 ) 



and 



(44) n = 



M* (1 + »jai)M*i + (1 + 77a 2 )M» 2 



Mbh = (j4 Ki + M2 K 2 ) 1,p + fV (aiM.i + tt2M. 2 ) 
M* {l+r]ai)M* 1 + (l + r)a 2 )M* 2 



respectively. In order to describe the effects on ry and ay of the radiative energy losses 
associated with gas dissipation, a fraction (1 + f)t] of the gas internal energy U g of each 
progenitor is subtracted from the total energy budget of the merger-product, consistent 
with the previous assumptions. Thus, the final total energy of the remnant is 

(45) E = E 1 + E 2 - 7] (1 + /) (aiJf.i + a 2 K* 2 ), 

and from the identities in eq. (37) one easily obtain 



(46) 
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1 _ {M*i + MgA 2 Ax fM* 2 + M g2 \ 2 A 2 



and 
(47) 

rv \M* + M gas J r V i ' \M* + M gas J r V2 ' 
where 

(48) Al = i + £+P*!l, 

and a similar expression holds for A 2 . In a dry (n = 0) merging Ai = A<i = 1, so that 

(49) mm a vl ,CT V2 < a v = — — < max <r vl , a V2 , 

(1 + ai)M*i + (1 + a2jM* 2 

i.e., the virial velocity dispersion of the merger-product cannot be larger than the maxi- 
mum velocity dispersion of the progenitors. Instead, A > 1 in the case of "wet" (r/ > 0) 
merging, and the resulting cry is larger than in the dry case, possibly larger than the 
maximum velocity dispersion of the progenitors. A similar argument shows that in the 
presence of gas dissipation the new ry increases less than in the dry case. Note that the 
above conclusions are obtained under the hypothesis of parabolic merging. If mergers 
involve galaxies on bound orbits, the additional negative energy term in equation (45) 
would lead to an increase of ay also in equal-mass dry mergers. The analysis of this case, 
and the question of how much fine-tuned the properties of the progenitor galaxies should 
be with their binding orbital energy in order to reproduce the SLs, are not discussed 
further in this review (for additional discussion, see, e.g., [241, 242]). 

The previous formulae are now applied to study repeated parabolic merging on a pop- 
ulation of Es. The merging spheroids are extracted by means of Monte-Carlo simulations 
from different samples of seed galaxies, the parent galaxies removed from the popula- 
tion and replaced by their end-product, whose properties arc determined by using the 
relations given above. The initial population of seed galaxies is obtained by random ex- 
traction of the stellar mass M* from the SDSS luminosity function in the Infrared ([265]), 
under the assumption of constant stellar mass-to- light ratio Y». For each galaxy mass, 
the corresponding central velocity dispersion o is fixed according to the FJ, and the 
effective radius R c is assigned from the FP ([74, 79]), while the observational properties 
of the resulting systems are determined from eq. (39). 

In the first scheme, the seed galaxies span only a narrow mass range (a factor of 
~ 5): in this case we then study whether massive ellipticals and the observed scaling 
relations can be produced by repeated mergers of low-mass spheroidal systems. In the 
second scheme the seed Es follow the observed scaling relations over their whole observed 
mass range (<~ 10 3 ), and so one explores whether repeated merging events preserve or 
destroy these relations. Figure 19 shows the results in the case of dry parabolic merging. 
The mass interval spanned by the progenitors is indicated by the two vertical ticks, the 
end-product positions are represented by dots, and the observed SLs are represented by 
dotted lines. It results that massive Es cannot be formed by parabolic dry mergers of low- 
mass spheroids only, because they would be characterized by exceedingly large values of 
R c and almost mass independent values of a , in agreement with the results of [31, 236]. 
Mergers with gas dissipation (not shown here) produce more realistic objects and the 
observed SLs are satisfied (even though with large scatter) by the new galaxies, up to a 
mass increase of a factor of 10 2 with respect to the smallest seed galaxies. However, new 
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Fig. 19. - Synthetic scaling relations produced by parabolic dry mergers. Seed galaxies span 
a limited mass range (indicated by the heavy vertical ticks) and random re-merging events 
are repeated until a factor 10 3 increase in mass is reached (see text for details). Dotted lines 
represent the observed scaling relations. All quantities are normalized to the properties of the 
lowest mass seed galaxy (from [251]). 



galaxies characterized by a mass increase factor >10 2 are mainly formed by mergers of 
gas poor galaxies that already experienced several mergers, and so they strongly deviate 
from the observed SLs. 

This first exploration therefore reveals that parabolic merging of low mass galaxies 
only is unable to produce Es obeying the observed scaling laws, even when allowing 
for structural weak homology in a way consistent with the edge-on FP. However, gas 
dissipation plays an important role in gas rich merging and, remarkably, the resulting 
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Es appear to be distributed as the observed SLs, as far as sufficient amounts of gas 
are available. Quite obviously, the problem of the compatibility of the properties of 
such merger-products with other key observations, such as the color-magnitude and the 
metallicity- velocity dispersion relations, and the increasing age of the spheroids with their 
mass (e.g. [58, 266]) remain opens. In the second scheme, the masses of the seed galaxies 
span the full range covered by ordinary Es (~ 10 3 ) and their characteristic size and 
velocity dispersion follow the observed SLs: in practice, we study the effect of merging 
(dry and wet) on already established SLs. The main result is that now, at variance with 
the results of the first scenario, the SLs remain almost unaffected by the merging, both 
in their slope and scatter. In particular, the Mbh-Co relation (with p = 2) is preserved, 
even though we are in a dry merging regime. The only detectable deviations from the 
observed SLs, for the same reasons discussed above, are found for Es with masses larger 
than the most massive galaxies in the original sample (marked by the two vertical ticks 
in Fig. 20). 

Why do mergers preserve so well the scaling relations? The reason is simple: by 
construction in a population of galaxies spanning the whole mass range observed today 
and distributed according to the observed SLs, mergers in general involve a "regular" Es, 
with values of R e and a as observed. These mergings act as a "thermostat" , maintaining 
values of R c in the observed range and increasing the virial velocity dispersion, thus 
contributing to preserve the SLs at increasing mass. Only when the produced galaxies 
are so massive that no regular galaxies of comparable mass are available, the new merger 
products deviate more and more from the SLs. This behavior becomes extreme in the case 
of repeated mergers in a galaxy population spanning a restricted mass range. Therefore, 
while Es cannot be produced by the merging of low mass spheroids only (as already 
pointed out by e.g. [31, 236, 240]) the observed SLs, once established by some other 
mechanism, are robust against merging. 

4 2. Cooling flows and SMBH feedback. - Quite often, when discussing the origin of 
the SLs between SMBHs and their host galaxies, the attention focuses on the galaxy 
formation mechanism, while almost no attention is payed to the subsequent several Gyrs 
of evolution (with the exception of possible merging events). This is probably due, 
quite surprisingly, to the common misconception that Es, once formed, are just "dead 
and boring red objects". However, as is well known for the past 30 years by stellar 
evolutionists and by the "cooling flow" community, this is just wrong (see Sect. 2). In 
fact, the mass return rate from the passively evolving stellar population sums up to a 
sizable (>10 — 20%) fraction of the galaxy stellar mass and is the main ingredient of the 
cooling flow model (and its variants). Clearly, some very efficient feedback mechanism 
supressing the cooling flow must be present in Es otherwise, in addition to the problems 
listed in Sect. 2, the SMBHs would have masses a factor ~ 100 greater than those 
observed and a luminous QSO should be present at the center of all Es. 

A possible (partial) solution to the cooling flow (and indirectly to the SMBHs growth) 
problems was proposed in [151, 267, 268] by considering the effect of SNIa heating of the 
galactic gas, and exploring the time evolution of gas flows by using hydrodynamical 
numerical simulations. It was found that while SNIa input sufficed for low and medium- 
luminosity Es to produce fast galactic winds, the inner parts of more massive spheroids 
would nevertheless host inflow solutions similar to cooling flows. This is because, while 
the number of SNIa per unit optical luminosity is expected to be roughly constant in 
ellipticals, the gas binding energy per unit mass increases with galaxy luminosity, as 
dictated by the FJ relation. However, as already recognized by [151], the mass budget 
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Fig. 20. - Synthetic scaling relations for the merger-products of up to 6 dry major mergers 
of galaxies extracted from a population that follows the observed scaling laws. Lines are as 
in Fig. 19 and all quantities are normalized to the properties of the lowest-mass seed galaxy 
(from [251]). 
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problem would still affect medium-large galaxies, putative hosts of luminous cooling 
flows. Thus, a concentrated feedback source is a very promising solution for a variety 
of problems, and the central SMBH is the natural candidate, by its mass and by its 
location, through a combination of mechanical and radiative feedback mechanisms. Some 
calculations have allowed for a physically motivated AGN feedback (e.g. [269]- [275]), 
and the computed solutions are characterized by relaxation oscillations. Energy output 
(radiative or mechanical) from the central SMBH pushes matter out, the accretion rate 
drops and the expanding matter drives shocks into the galactic gas. Then the resulting 
hot bubble ultimately cools radiatively and the resulting infall leads to renewed accretion. 
The cycle repeats. Among the computed models that studied the interaction between 
AGN feedback and galactic cooling flows, those of [270, 271, 275] focused on the effects 
of radiative heating on galactic gas flows. In fact, if one allows the radiation emitted from 
the accreting SMBH to interact with and heat the galactic gas, one solves the cooling 
flow problem in Es, and the feedback produces systems that are variable but typically 
look like normal ellipticals containing hot gas. They sometimes look like incipient cooling 
flows and rarely, but importantly, appear like quasars. Interestingly, observations seem 
to support this scenario (e.g. [276]). 

In [270, 271], however, a major uncertainty remained about the typical QSO spec- 
trum to adopt, in particular the high energy component of that spectrum, which is most 
important for heating the ambient gas. Thus, a simple broken power law was adopted for 
the spectrum with a range of possible values of the Compton temperature - from 10 72 K 
to 10 9 ' 5 K - with most of the emphasis of the paper being on the higher temperatures. 
Subsequent work by Sazonov, Ostriker & Sunyaev ([277], see also [279]), which assessed 
the full range of observational data of AGNs and computed their Spectral Energy Dis- 
tribution, concluded that the typical equilibrium radiation temperature was narrowly 
bounded to values near 10 7 3 K, i.e., of the order of 2kcV. This value is still above the 
virial temperature of all galaxies and, most importantly, well above the central tempera- 
ture of the cooling flow gas. As noted in [22], there is a rather large compensating effect 
also not included in [271]: gas heated by radiation with a characteristic temperature near 
10 7 K is heated far more effectively by absorption in the atomic lines of the abundant 
metal species than by the Compton process. In particular, [22] provide a fitting formula 
for the Compton plus photoionization and line heating/cooling that was implemented 
into the numerical code developed by Ciotti & Ostriker. 

Here I briefly describe the main results of [275]. Consistent with HST observations, 
which have shown that the central surface brightness profile is described by a power- 
law (see Sect. 2.1.3) as far in as it can be observed (i.e. to ~ 10 pc for Virgo ellip- 
ticals), the stellar component is described by a Hernquist model. In addition, optical 
(e.g. [214], [280]-[282]), and X-ray (e.g. [283, 284]) based studies typically find that lumi- 
nous matter dominates the mass distribution inside the effective radius R ej while dark 
matter begins to be dynamically important at 2 — 3i? c , with common values of the to- 
tal dark-to-luminous mass ratio 1Z = Mh/M* in the range 1<1Z<6. Finally, theoretical 
(e.g. [206, 187, 209, 285] and observational (e.g. [286]) arguments support the idea that 
similarly to the stellar profiles, also the radial density distribution of the dark halos is 
described by a cuspy profile. Following these empirical indications, in the numerical sim- 
ulations of [275] the stellar and dark matter distributions are taken to be of the form of 
eq. (28), with Mh = and — [3r* the halo total mass and scale- length; dynamical 

and phase-space properties of the resulting two-component Hernquist models are given 
in [191]. The physical scales of the model are then fixed so that the model satisfies the 
FJ and the edge-on FP. Once the density profile is fixed, the velocity dispersion profile is 
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computed by solving the Jeans equations, as this quantity is an important ingredient in 
the energy budget of the gas flows, namely the thermalization energy of the stellar mass 
losses. The stellar mass losses are computed following the detailed prescriptions of stellar 
evolution, while the SNIa rate (the time dependence of which is unfortunately still quite 
uncertain) is normalized to recent empirical estimates ([287, 288]). Each SNIa event re- 
leases Esn — 10 51 erg of energy and 1.4M Q of material in the Interstellar Medium. The 
simulations allow for star formation, which cannot be avoided when cool gas accumulates 
in the central regions of Es. In the new population, described by a Salpeter IMF, the 
associated total number of Type II Supernovae and the stellar mass return are also com- 
puted, together with the radial distribution of the optical and UV luminosity per unit 
volume of the new stars. The accretion disk around the SMBHs (feeded by hydrodynam- 
ical processes) is incorporated as a set of ordinary differential equations describing the 
istantaneous mass of gas and stars (low and high mass, and remnants) in the disk, as well 
as the disk optical and UV luminosity. Finally the mass accretion on the central SMBH 
in terms of the so-called a-viscosity prescription is obtained, leading to the bolometric 
accretion luminosity 

(50) L BH (t) = eM BH (t)c 2 . 

For the radiative efficiency the standard value e = 0.1 is assumed (as suggested by ob- 
servations, e.g. [36, 37, 289]), but a generalization to include an ADAF-like efficiency 
([271, 290]) is also explored. The radiative heating and cooling produced by the accre- 
tion luminosity are computed from the formulae provided in [22] , which describe the net 
heating/cooling rate per unit volume E of a cosmic plasma in photoionization equilibrium 
with a radiation field characterized by the average quasar Spectral Energy Distribution 
given in [277], for which the associated spectral (Compton) temperature is Tx — 2 keV. 
In particular, Compton heating and cooling, bremsstrahlung losses, line and recombi- 
nation continuum heating and cooling are taken into account. Radiation pressure on 
the Interstellar Medium (via electron scattering, dust opacity, photoionization opacity) 
produced by accretion luminosity and by stellar light (consequence of star formation) is 
obtained numerically by solving the two lowest spherically symmetric moment equations 
of radiative transfer in the Eddington approximation (e.g., [291]). 

The evolution of the galactic gas flows is obtained by integrating the time-dependent 
Eulerian equations of hydrodynamics, 

(51) ^7 + V-(pv) = ap*+p„-p+ 



(52) — + V • (mv) = -( 7 - 1)VE - Vp rad + pg - m+, 



dE 

(53) — + V- (Ev) = -(7- l)EV- w + H-C 

„2 i Q„2 



^iM^+s, + E ll -Ei. 



where p, m, and E are the gas mass, momentum and internal energy per unit volume, 
respectively, and v is the gas velocity (for a complete description of the source and sink 
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terms, for a discussion of the adopted assumptions, and for an outline of the numerical 
integration scheme, see [275, 292]). The initial conditions are represented by a very low 
density gas at the local virial temperature. The establishment of a high-temperature gas 
phase at early cosmological times is believed to be due to a "phase-transition" when, as 
a consequence of star formation, the gas-to-stars mass ratio was of the order of 10% and 
the combined effect of SNIa explosions and AGN feedback became effective in heating 
the gas and driving galactic winds. Several theoretical arguments and much empirical 
evidence, such as galaxy evolutionary models and the metal content of the Intracluster 
Medium, support this scenario (e.g. [24, 293]). 

What are the main phases of the model evolution? As an illustration, Figs. (21)-(25) 
refer to a model with an initial stellar mass M* = 4.6 x 10 11 M ( t ) , effective radius R e = 6.9 
kpc and a — 260 km s _1 (chosen to lie on the FP), total dark-to- visible mass ratio TZ = 5 
and dark-to-visiblc scale-length ratio j3 — 5.22 (corresponding to an identical amount of 
stellar and dark matter within the half-light radius). The initial SMBH mass follows the 
present day Magorrian relation, i.e., Mbh — 10 _3 M*. After a first evolutionary phase in 
which a galactic wind is sustained by the combined heating of SNIa and thermalization of 
stellar velocity dispersion, the central "cooling catastrophe" begins. In the absence of the 
central SMBH a "mini-inflow" would be then established, with the flow stagnation radius 
(i.e., the radius at which the flow velocity is zero) of the order of a few hundred pc to a few 
kpc. These decoupled flows are a specific feature of cuspy galaxy models with moderate 
SNIa heating ([268]). However, after the central cooling catastrophe, the feedback caused 
by photoionization and Compton heating strongly affects the subsequent evolution, as 
can be seen in Fig. 21. The bolometric luminosity (top panel) ranges between roughly 
0.001 to 0.1 of the Eddington limit at peaks in the SMBH output but, since obscuration 
is often significant, the optical accretion luminosity as seen from infinity can be much 
lower (bottom panel). The major AGN outbursts are separated by increasing intervals of 
time (set by the cooling time) and present a characteristic temporal substructure, whose 
origin is due to the cooperative effect of direct and reflected shock waves. At t ~ 8 Gyrs 
the SNIa heating, also sustained by a last strong AGN burst, becomes dominant, a global 
galactic wind takes place, and the nuclear accretion switches to the optically thin regime. 

Remarkably, the coronal X-ray luminosity Lx, due to the hot galactic atmosphere, 
falls in the range commonly observed in massive Es, with mean values lower than the 
expected luminosity for a standard cooling flow model. It is also found that a large 
fraction of the starburst luminosity output occurs during phases when shrouding by 
dust is significant. An important quantity associated with the time evolution of the 
various luminosities is their duty cycle (for the operational definition see [271, 275]), 
whose temporal evolution is showed in Fig. 22: of course, the duty-cycles of starburst 
optical and UV luminosities ar larger and less fluctuating than those of the AGN, and 
overall they are in agreement with the observations (e.g. [294]). Of particular interest 
for the present discussion is the evolution of the mass budget of the model. In Fig. 23 
the time evolution of some of the relevant mass budgets of the model, both as time- 
integrated properties and instantaneous rates, is shown. At the end of the simulation 
the total Interstellar Medium mass in the galaxy is ~ 5 x 10 8 M Q , while the SMBH 
mass reaches a final mass of ~ 7 x 10 8 M Q . a model with a smaller initial SMBH mass 
would accrete less, thus maintaining the Magorrian relation even better. The SMBH 
mass accretion rate strongly oscillates as a consequence of radiative feedback, with peaks 
of the order of 10 or (more) M Q /yr, while during the final, hot-accretion phase the 
almost stationary accretion rate is <10~ 4 M Q /yr: this value is close to the estimates 
obtained for the nuclei of nearby galaxies ([165]). Note that in the last 6 Gyrs the SMBH 
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Fig. 21. - Dotted lines are the bolometric accretion luminosity (top panel) and the optical 
SMBH luminosity corrected for absorption (bottom panel). The almost horizontal solid line is 
the Eddington luminosity L^dd (from [275]). 
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Fig. 22. - Time evolution of duty cycles. Top panel: duty-cycle of I/bh,uv (solid) and I/g H opt 
(dotted); the top axis shows the corresponding redshift. We see that these systems would be 
observed from afar in the (rest-frame) optical or UV as quasars several percent of the time. 
Bottom panel: duty-cycle of the starburst (solid), L^pt (dotted), of the Interstellar Medium 
X-ray luminosity (computed in a volume excluding the inner 100 pc), and of the recycled IR 
luminosity Lir (from [275]). 
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virtually stops its growth, while the Interstellar Medium mass first increases, due to the 
high mass return rate of the evolving stellar population, and then decreases, due to the 
global galactic wind induced by SNIa. During the entire model evolution, more than 
10 10 - 5 M Q of recycled gas are added to the Interstellar Medium from stellar mass losses. 
Approximately 2.1 x 1O 1O M have been expelled as a galactic wind, while ~ 1.4 x 10 10 M Q 
are transformed into new stars, so that only 0.7% of the recycled gas is accreted onto the 
central SMBH. The central paradox of the mass budget is thus resolved. An identical 
model without SMBH feedback, but with the same star formation treatment of the model 
described above, would produce a SMBH of final mass >1O 1O M0, while the total mass 
in new stars would be reduced to ~ 3 x 10 9 M Q . The star formation rate during the 
periods of feedback dominated accretion oscillates from 0.1 up to several hundreds (with 
peaks near 10 3 ) M Q yr" 1 , while it drops monotonically from 10 _1 to <10~ 3 M Q yr^ 1 
in the last 6 Gyrs of quiescent accretion. These violent star formation episodes (with 
SMBH accretion to star formation mass ratios <~ 10~ 2 or less) are induced by accretion 
feedback and are spatially limited to the central 10 — 100 pc; thus, the bulk of gas flowing 
to the center is consumed in the star burst. These findings arc nicely supported by recent 
observations (e.g [295, 296]). Note that the "age" effect of the new stars on the global 
stellar population of the galaxy is small, as the new mass is only 3% of the original stellar 
mass and it is virtually accumulated during the first Gyrs (see Fig. 23). 

In Fig. 24 the final spatial density profile of the system is shown, together with its 
projection and the best-fit obtained with the Sersic law. As expected, the profiles show 
an increase of the best-fit Sersic parameter to, due to the mass accumulation in the 
central regions. Remarkably, the final value of m is within the range of values commonly 
observed in ellipticals: however, in the final model we note the presence of a central (~ 30 
pc) nucleus originated by star formation which stays above the best fit profile, similar to 
the light spikes characterizing "nucleated" galaxies. 

Figure 25 shows the temperature and density in the central regions of the model: note 
how the SMBH bursts heat the central gas, causing the density to drop, and launching 
gas at positive velocities of the order of thousands km s _1 . The Compton temperature Tx 
is the horizontal dashed line; during the bursts, the local gas is heated above this limit. 
As was already found in [271], the galaxy cooling catastrophe starts with the formation 
of a cold shell placed around the galaxy core radius: however, in the present models (see 
also [268]), the cooling catastrophe happens at significantly earlier times, because of the 
higher central stellar density and of the different time dependence and amount of SNIa 
explosions. This cycle of shell formation, central burst, and expanding phase, repeats 
during all the bursting evolution, along the lines described in detail in [271]. Finally, 
when the SNIa heating per unit gas mass becomes dominant over the decline of fresh 
mass input from evolving stars, the galaxy hosts a wind, and the accretion becomes 
stationary without oscillations; the central SMBH radiates at <~ 10 _5 LEdd (e.g. [297]). 

4'3. Summary. - In this Section I briefly addressed an astrophysical question possibly 
as important as the origin of the galaxy and SMBH Scaling Laws, i.e. the robustness of the 
relvant SLs against physical phenomena (in principle) able to destroy them. In particular, 
I focused on galaxy merging and SMBH accretion at the center of cooling flows fed by 
stellar mass losses, which represents the major contributor to the Interstellar Medium in 
the life of Es after their formation. The main conclusions can be summarized as follows: 

1) Parabolic dry merging in a population of low mass spheroids leads to massive Es 
that fail the FJ and Kormcndy relations, being characterized by low velocity dispersion 
and very large effective radii. Parabolic wet merging in the same population of low mass 
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Fig. 23. - Mass budget evolution. Top left panel: total hot gas mass in the galaxy (within 10R a , 
Mism, dotted line), and accreted mass on the central SMBH (AA/bh). Top right: mass lost as 
a galactic wind at 10i? c (AA/ gw , dashed line), and total mass of new stars (AM*). Bottom left 
panel: global mass return rate from the evolving stellar population and mass accretion rate on 
the central SMBH (Mbh)- Bottom right: galactic wind mass loss rate at 10J? e (Af gw , dashed 
line) and instantaneous, volume integrated, star formation rate. Note that for t > 8 Gyrs 
the mass lost as a galactic wind is almost coincident with the mass input from evolving stars 
(from [275]). 
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Fig. 24. - Dotted lines are the projected surface density of the model shortly after the beginning 
of the simulation (t = 2.5 Gyr, z ~ 2.6, bottom panel) and at t = 13.5 Gyr (z — 0, top panel), 
normalized to the surface density at the effective radius. Solid lines are the best-fit Sersic law. 
The effective radius contracts from ~ 9.2 kpc to ~ 8.4 kpc, while the surface density E c increases 
from ~ 3 x 10 22 to ~ 3.6 x 10 22 protons per cm -2 (from [275]). 



progenitors leads to galaxies in better agreement with the observed scaling relations, as 
long as enough gas for dissipation is available. 

2) The edge-on structure of the FP is surprisingly preserved (except in the case of 
numerous minor mergers with galaxies on radial orbits). In fact, in the last case the devi- 
ations from the FJ and Kormendy relations do not compensate. Another manifestation 
of the problem encountered by head-on minor mergers is that the Sersic parameter m 
describing the mass profile of the end-products decreases for increasing mass, at variance 
with observations. In all cases, galaxies remain in the populated region of the face-on 
FP. Incidentally, points 1) and 2) show that the FJ and Kormendy relations, despite 
their larger scatter, are stronger tests for merging than the edge-on FP, thus providing a 
powerful way to investigate the assembly history of massive elliptical galaxies. 

3) Parabolic dry or wet mergers in a population of galaxies following the observed 
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Fig. 25. - Top panel: gas velocity at 5 pc from the SMBH. Note how the SMBH growth affects 
the lower envelope of velocity values. Bottom panel: Gas number density (dotted line, scale on 
the right axis) and temperature at 5 pc from the SMBH (solid line, scale on the left axis). Low- 
temperature, high-density phases end when accretion luminosity Lbh increases sharply heating 
the ambient gas to a high-temperature, low-density state. The horizontal dashed line is the 
model Compton temperature Tx = 2.5 x 10 7 K (from [275]). 
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scaling laws over the full mass range populated today by stellar spheroids, preserve the 
Kormendy, FJ, and edge-on FP remarkably well. The reason of this behavior is due 
to the presence, in the merger population, of galaxies with velocity dispersion increasing 
with galaxy mass. Thus, massive Es cannot be formed by (parabolic) merging of low 
mass spheroidal galaxies, even in presence of substantial gas dissipation. However, the 
observed scaling laws of Es, once established by galaxy formation, are robust against 
merging. 

4) Under the reasonable hypothesis that the derived values of a arc not strongly 
affected by the dynamical evolution of binary BHs, dissipationlcss merging, while in ac- 
cordance with the Magorrian relation, fails to reproduce the Mbh-Co relation. Curiously, 
by allowing for substantial emission of gravitational waves during the BHs coalescence, 
the Mbh-Co relation is reproduced, but the Magorrian relation is not. In the case of gas 
dissipation, the resulting Mbh-^o relation is in better agreement with the observations, 
also in the case of significant mass loss (via gravitational waves) of the coalescing BHs. 
However, significant deviations from the observed scaling laws are expected for massive 
galaxies. 

5) The recycled gas from dying stars is an important source of fuel for the central 
SMBH, even in the absence of external phenomena such as galaxy merging, which are 
often advocated as the way to induce QSO activity. Radiative feedback from a central 
SMBH has dramatic effects on its mass growth: much of the recycled gas falling towards 
the galaxy center during the accretion events is consumed in central starbursts with a 
small fraction (of the order of 1% or less) accreted onto the central SMBH. Thus, the 
central starburst regulates the amount of gas available to be accreted onto the central 
SMBH. If we did not allow for the (AGN feedback induced) central star formation, the 
SMBH would grow to be far more massive than seen in real galaxies. 

5. Toward a unified picture? 

After discussing how the SLs can be maintained over cosmological times, we have 
now to address plausible scenarios aimed at explaining how these laws are established. 
In this Section I discuss three seemingly unrelated arguments, that will be connected in 
Sect. 6. In particular, in Sect. 5.1 I discuss the growth of SMBHs, in Sect. 5.2 I describe 
the results of numerical simulations of fast (dissipationlcss) collapse in pre-existing dark 
matter halos, and finally in Sect. 5.3 I consider the remarkable SLs inprinted in dark 
matter halos produced by cosmological simulations. 

5T. SMBH and spheroid growth. - It is now universally accepted that the SLs involv- 
ing the central SMBHs and the global properties of the host spheroids are established at 
the epoch of galaxy formation. In other words, the Magorrian relation reveals that the 
bulk of SMBH fueling in AGNs must be associated with star formation in the spheroidal 
components of their host galaxies (e.g. [18, 29, 30, 37], [316]-[318]). This argument is 
rich in consequences, as it naturally links galaxy formation to the SMBH growth. A first 
quantification of this link is dprovided by the so-called Soltan argument ([289]), which in 
its more recent applications (e.g. [36, 37]) shows that the integrated luminosity emitted 
by QSOs over the life of the Universe nicely matches the present day total mass of SMBHs 
at the center of Es, for efficiencies of <~ 0.1. In particular, the approach in [37] is based on 
two working hypotheses, i.e. that 1) spheroid star formation and BH fueling are - at any 
time and in any system - proportional to one another with the proportionality constant 
independent of time and place, and that 2) the SMBH accretion luminosity always stays 
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near the Eddington limit when the QSO is in the luminous phase and the BH does not 
produce any radiation in the "oft" state (e.g. because accretion is suppressed). Assump- 
tions 1) and 2) are then coupled with three observational inputs, namely the present-day 
luminosity function of spheroids, where the number of spheroids per unit volume with 
rest-frame B-band luminosities in the interval (Lb, Lb + oLLb) is given by $s^b, with 

3>S« I Lb \ ( L\ 



(54) W = E^x f 1 exp 



i=i 



Ls*i \Ls*i J \ ^S*i 



and the different indices correspond to the different galaxy types mentioned in Footnote 
1. The second ingredient is the luminosity function of QSO 

/rr\ ^ IT \ ®Q*/Ln*(z) 

(55) &q(Lq,z) - 



[L Q /L Q *{z)]0i + [L Q /L Q *(z)]^ 
where the characteristic luminosity Lq* in the rest-frame B band is 

.! e^(l + e« z *) 



(56) L Q »(z) = L Q ,(0)(l + z) a «- 



(e.g. [33, 34]). Finally, the third ingredient is the Magorrian relation, obtained by com- 
bining the FJ with the Mbh — c relation 

(57) 0.016^. 

(see also Sect. 2). Note that under assumption 2) it is expected that the redshift evolution 
of the QSO emissivity and of the star formation history in spheroids should be roughly 
parallel to each other: indeed, numerical simulations of feedback-modulated accretion 
flows (radiative, as in [271, 275], or mechanical, as in [269, 319, 320]) show that the 
accretion luminosity during short episodes of bursts stays near the Eddington value (see 
Fig. 22). 

For illustrative purposes, let us first consider a population of A g identical galaxies 
over the Hubble time tn, each of which today (i.e. at t = tn) hosts a spheroid of mass 
Ms, and a SMBH of mass Mbh- Let us further assume that during the entire time 
elapsed from to £h, each SMBH had only two states: it was either "on" or "off". We 
identify the "on" state as the active quasar phase, and we define the duty cycle /q as the 
fraction of the time each BH spends in the "on" state. At any given time, the number of 
active quasars is then A^q = fc}N g . In the "on" state, the SMBH grows by accretion at 
the rate Mbh, and shines at the (bolometric) luminosity Lq with radiative efficiency e, 
defined as the fraction of the rest mass energy of the infalling gas converted to radiation. 
The remaining fraction (1 — e) of the rest mass then leads to the growth of the BH mass 
([36]). Simple algebra shows that 

( 58 ) _J_ E Q fQt H L Q N g _ t H L Q N Q 



MT hC 2 M BH A k c 2 Mbh^c 2 



Here c is the speed of light; the numerator represents the total light emitted by all 
SMBHs, and the denominator represents the total mass in SMBHs today. In the third 
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equality, we have used Nq — fqN g . Note that the last term involves only quantities 
that are, in principle, directly observable and that it is independent of the duty cycle. 
Equation (58) describes the entire galaxy population, but a similar equation applies to 
individual galaxies: Lq/q^ — eMbhc 2 . This last expression does have a dependence on 
the duty cycle, which can therefore be written as 

(59) M = £ = ^ 



The above argument demonstrates that the radiative efficiency can be obtained indepen- 
dently of the duty cycle and that the duty cycle can be obtained in two different ways, 
based either on the number or on the characteristic SMBH mass of quasars. By using 
bona-fide observed values for the different quantities involved in the above equations, 
in [37] it is found that e = 0.071, in good agreement with the result of [36]. 

The duty cycle defined in eq. (59) can be generalized to a population of evolving 
galaxies under the assumption that the duty cycle does not vary with time but it is a 
function of luminosity. This is done by defining the average duty cycle of all quasars 
above luminosity Lq(x, 0), 

where Lq(x, 0) is such that QSOs at redshift z = brighter than Lq(x, 0) emit a fraction 
x of the total quasar light Lq = J °° dL L$q(L, 0). Likewise, if Mbh(x) is such that all 
SMBHs more massive than Mbh(x) sum up to the same fraction x of the total SMBH 
mass at z = 0, then the last term in eq. (59) becomes 

(61) \Jq,m)x 



J* °° L Q {x,t)dt 

It is found that (/q,n)o.i — 0.008 and (/q ; n)o.9 — 0.05, and that the two methods agree 
well on the high mass end, while (/q,m) is systematically lower by a factor of ~two 
towards the low-mass end. These results for the duty cycle are in good agreement with 
the values reported in [36], [321]- [323], and with theoretical expectations (e.g [270, 271, 
275]). 

In a complementary aproach to the previous phenomenological investigation, it is nat- 
ural to attempt the modeling of the simultaneous growth of Es and their central SMBHs. 
While numerical simulations of galaxy formation are becoming richer and richer in the 
input physics, simple "toy models" are still helpful, because of their specific capability 
to incorporate in a intuitive way the core physics of the investigated process. As an 
example, here I describe the results obtained with the model explored in [22, 272] (see 
also [15, 17], [19]- [21], [30, 250, 321, 318]). The basic idea of [22] is that the UV and 
high energy radiation from a typical quasar can photoionize and heat a low density gas 
up to an equilibrium Compton temperature (T c w 2 x 10 7 K) that exceeds the virial 
temperatures of giant ellipticals. Note that the radiative output is not the only, nor 
even necessarily the dominant mechanism whereby feedback from accretion onto central 
SMBHs can heat gas in Es. For example, [269, 319] have stressed that the mechanical 
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input from radio jets will also provide a significant source of energy, and much detailed 
work has been performed to follow up this suggestion. 

An empirical indication that radiative feedback from the accreting SMBHs can in- 
deed lead to final masses following the observed scaling with the host system velocity 
dispersion has been obtained by Sazonov et al. ([22, 277]). In fact, Sazonov et al. ([277]) 
computed for the observed average quasar Spectral Energy Distribution the equilibrium 
temperature T oq of gas of cosmic chemical composition as a function of the ionization 
parameter £ = Lbh/(w 2 ), showing that 

( 10 4 K for £ < 100: 

(62) T oq (0« { 2xl0 2 £K; 

[2xl0 7 K for £>5xl0 4 , 

where n is the hydrogen number density and r is the distance from the SMBH. In practice, 
T eq is the temperature at which heating through Compton scattering and photoionization 
balances Compton cooling and cooling due to continuum and line emission. Suppose now 
that the gravitational potential experienced by the gas is due to the central SMBH alone. 
Then the condition 

(63) 5 -kT cq (0 - ^JI^p > o 

(fj, is the mean molecular weight) roughly defines a situation where gas of density n, 
located at r, will be heated to T cq by the central radiation and blown out of the SMBH 
potential. Thus, for given Mbh, ^BH/^Edd and r, gas with density below a certain 
critical value, cannot accrete onto the SMBH. In terms of Lbh, this means that Bondi 
accretion ([278]) of gas at temperature T can be disrupted if Lbh is sufficiently high that 
T cq (R-B) ^ T, or equivalently L B h > L ciit , where 

(64) Lerft (T) = Z(T)Rl(T)n(R B ), 

£(T) is the ionization parameter corresponding to T oq = T, and finally 

, fi ,, r> GMBH^m p 1 M BH ( T 

(65) = 2jkT = 16P %10^ [Wk 

is the Bondi radius. 

The previous argument suggests that, before the SMBH grows to a certain critical 
mass M B H,crit, its radiation will be unable to heat efficiently the ambient gas, and accre- 
tion onto the SMBH will proceed at a high rate. Once the SMBH has grown to M B H,crit, 
its radiation will heat and expel a substantial amount of gas from the central regions of 
the galaxy. Feeding of the SMBH will then become self-regulated on the cooling time 
scale of the low density gas. Subsequent quasar activity will be characterized by a very 
small duty cycle (<~0.001), and the SMBH growth will be essentially terminated. On a 
more quantitative level, suppose that the galaxy density distribution is that of a singular 
isothermal sphere, with the gas density following the total density profile: 



50 



L. CIOTTI 



where M gas and M* are the gas mass and the stellar mass within some reference radius. 
Then, simple algebra shows that radiation from the central SMBH can heat the ambient 
gas up to the temperature 

(67) T cq « 6.5 x 103 I^L^J^L. f^s^V R 

V ' q ^Edd M gas 10 8 M Q \ a J 

The transition from rapid SMBH growth to slow, feedback limited SMBH growth is 
expected to meet the critical condition T cq = i] csc T* (where r^sc ^ 1 and T» is the 
temperature associated with the galaxy velocity dispersion a), so that 



(68) Mi 



Therefore, for fixed values of T] esc , LBH/^Edd and M gas /M* one expects MBH,crit oc <r 4 . 
It follows that the observed Mbh-Co relation will be established if 

(69) M&l ~ 3 x 10 ' 3 l bh 

T] csc Z/Edd 

To satisfy the observed Mbh-c relation, the gas-to-stars ratio is thus required to be 
relatively low and approximately constant for spheroids with different masses at the 
epoch when the SMBH reaches its critical mass. 

The approximately linear T cq (£) dependence is crucial to the above argument leading 
to the MBH.crit oc er 4 result. However, the T eq (£) function becomes strongly nonlinear 
outside the range 2 x 10 4 K < T cq < 10 7 K and a more general result can be obtained if 
we consider the exact curve T eq (£) from [277]. The situation is summarized in Fig. 26. 
It is perhaps interesting that the range of masses shown in Fig. 26 for which Mbh oc 
<t 4 is obtained from considerations of atomic physics (and the observed AGN spectra) 
corresponds closely with the range of masses for which this power law provides a good 
fit to the observations. 

Starting from the previous results it is possible to study the SMBH-galaxy co-evolution, 
by using a physically-motivated one-zone model. The model differential equations for the 
gas mass budget M gas of the galaxy adopted in [22] are 

(70) M gas = M inf — M* + M re c - Mbh - M esc , 

where the quantities on the r.h.s. describe the cosmological infall on the forming galaxy, 
the amount of gas subtracted by star formation, the gas produced by evolving stars, the 
gas accreted on the SMBH, and finally the gas lost as a galactic wind when the thermal 
energy of the Interstellar Medium is high enough to escape from the galaxy potential 
well, respectively. More in detail, 



(71) 



M inf = ^e-/-, M*= a ; M ^ -U, 

Tinf max(Tdyn,Tcool) 



rcci 



M rec = / M+{t')W*(t-t') dt', 
K JO 
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Fig. 26. - Thick solid line shows the predicted Mbh-cto relation resulting from the requirement 
that heating of the interstellar gas by radiation from the central SMBH at the Eddington limit 
be below the level required to drive the gas from the galaxy (T cq < T*, M gas /M = 0.003 and 
r/csc = 1). The thin solid line corresponds to M gaa /M = 0.0015 and ri esc — 2. The dashed line 
is the observed Mbh oc ai relation in the range 10 6 < Mbu/Mq < 10 9 . The dotted lines are 
M B h oc al and M B h oc al laws (from [22]). 



(72) M BH = M BH ,acc + Pbh,*M+, M BH ,acc = M inf (/ Edd M Edd , M B ), 

where is the stellar mass formation rate, and finally 



CSC 



(73) M c 
The dynamical time T dyn is defined as 



~ i * — '/esc-* * ) 
7"esc 

o, T< Vesc T*. 



(74) r dyn = — A 

where w c is some characteristic circular velocity of the dark matter halo. For an isothermal 
halo with 1-dimensional velocity dispersion a 

_ [im v a 2 _ nm p v 2 c 
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Moreover, the cooling time is evaluated as 

E 

(76) r cool = 7r , E c = n e n p A(T), 

Ec 

where A(T) is the gas cooling function, and E the gas internal energy per unit volume. 
The gas recycled by the evolving stellar population is obtained as a convolution between 
the instantaneous star formation rate and the kernel W*(i) derived from stellar evolu- 
tion ([151, 270, 271, 275]). In eq. (72) the Eddington and Bondi accretion rates are given 

by 

(77) M Edd = M B - 47ri?| PB c s , 
where 0.001 < e < 0.1, the Eddington luminosity is given by 

(78) L Edd = 1.3 xlO 46 ^^ ergs' 1 , 



and c s = y/kT/nnip is the (isothermal sound velocity). 

Finally, the thermal state of the gas is described by the rate at which (internal) energy 
per unit volume changes with time: 

(79) E = £?h,sn + £h,w + -Eh,agn — Ec + Eflf 



where the terms on the r.h.s. are the heating produced by supernova explosions, the 
thermalization heating of stellar winds, AGN heating, gas cooling and finally adiabatic 
heating and cooling due to galactic winds and cosmological accretion. 

The time evolution of the quantities shown in Fig. 27 refers to a galaxy model char- 
acterized by R c = 4kpc and a halo (constant) circular velocity of 400kms _1 ; the total 
mass of the gas infall is lO n M , and the characteristic infall time is 2 Gyr. The initial 
black hole mass is assumed to be 10M Q . The complete description of the toy-model 
behavior for different choices of the input parameters is given in [22]. Here I just recall 
that after an initial "cold" phase dominated by gas infall, as soon as the gas density 
becomes sufficiently low, and correspondingly the cooling time becomes longer than the 
dynamical time, the gas heating dominates, and the galaxy switches to a "hot" solution. 
The gas mass/stellar mass ratio (<~ 0.003) at that moment is remarkably near to the 
value inferred from the argument leading to the right Mbh-0o relation. Note also how 
the gas content of the present day galaxy model and the final black hole mass are in nice 
agreement with observations. 

An interesting experiment is obtained by reducing the circular halo velocity and the 
infall mass in the reference model. In these cases, galactic winds, powered by supernova 
heating, are favoured, i.e. small galaxies lose their gas content easily, in accordance with 
the predictions of hydrodynamical simulations and as expected from the Mg2-<7 and 
the FJ relations. Remarkably, the transition to the hot phase of these models happens 
for Mg as /M* ~ 0.01, similarly to the behavior seen in the case of the more massive 
spheroid in the reference model. This sort of cooperation between AGN feedback and 
stellar energy injection, i.e., the fact that in general substantial galactic winds are due to 
stellar heating, and are reinforced by the presence of the central A GN, was already found 
in numerical simulations. 
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Fig. 27. - Time evolution of the model gas temperature (solid line) . The model virial temperature 
is represented by the dotted line, while the dashed line represents the "escape" temperature (here 
assumed 2T„). Note the strongly fluctuating temperature in the initial "cold" phase (from [22]). 



An important and apparently robust conclusion that can be drawn from these sim- 
ulations is that inevitably stellar heating leads to a transition from cold to hot solution 
when the gas-to-star mass ratio drops to of order 1 per cent or somewhat less. Now, since 
a gas fraction of this order is required for the radiative feedback from the central SMBH 
to limit its growth at the mass obeying the observed Mbh-c relation, it is tempting 
to suggest that the SMBH reaches its critical mass, determined by radiative feedback, 
approximately at the epoch of transition from cold to hot galaxy phase. The near co- 
incidence of the gas fraction corresponding to the beginning of the hot galactic phase 
with that in eq. (69) required by our argument leading to the correct Mbh-c relation 
offers the possibility of the following evolutionary scenario. At the early stages of galaxy 
evolution, when the protogalactic gas is dense and cold, active star formation is accom- 
panied by the growth of a central SMBH. However, the black hole is not massive enough 
to produce a strong heating effect on the ambient, dense gas, even during episodes when 
it shines near the Eddington limit. This cold phase would be identified observationally 
with the Lyman Break Galaxies and bright submillimctcr galaxies, which are charac- 
terized by high star formation rate and moderate AGN activity. The cold phase ends 
when the gas-to-star mass ratio has been reduced to <~ 0.01, when the energy input from 
the evolving stellar population and possibly from the central SMBH heats the gas to 
a sub-virial temperature. The SMBH continues to grow actively during this transition 
epoch (that would be identified with the quasar epoch), because there are still sufficient 
supplies of gas for accretion, and soon reaches the critical mass (obeying the Mbh-Co 
relation), when the SMBH radiative output causes a major gas outflow. The subsequent 
evolution is passive and characterized by AGN activity with a duty cycle reduced by 
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a factor of ten to 0.001; this late phase would be identified with the present day Es, 
discussed in Sect. 4. 

Obviously, the results described above should not be ovcrintcrpretcd. As is common 
in studies based on a similar approach, the parameter space is huge (even though several 
input parameters are nicely constrained by theory and/or observations), so that the 
results of simulations of this kind should be interpreted more as indications of possible 
evolutionary histories than exact predictions. In particular, the toy-model cannot directly 
test the ability of radiative feedback to produce the right final SMBH mass. In fact, this 
can be done only using true hydrodynamical simulations. This is not surprising, because 
the toy-model, by construction, is a one zone model, and we already know that feedback 
mechanisms are strongly scale-dependent, in the sense that central galaxy regions react 
in a substantially different way with respect to the rest of the system. 

5"2. Collapse. - It is a well established fact that the end-products of dissipationless 
collapse reproduces several structural and dynamical properties of Es. For example, the 
pioneering work of van Albada ([60]) showed that the end-products of cold collapse have 
projected density profiles well described by the R 1 ^ 4 de Vaucouleurs law, radially decreas- 
ing line-of-sight velocity dispersion profiles, and radially increasing velocity anisotropy, 
in agreement with observations of Es ([325, 326]). More recently, dissipationless col- 
lapse has been studied in greater detail thanks to the advances in N-body simulations 
(e.g. [327]- [330]). These studies show that a smooth final density distribution with i? 1 / 4 
projected mass profile is produced when the initial conditions are cold, extended, and 
clumpy in phase-space. From the astrophysical point of view the dissipationless collapse 
([47]) was introduced to describe a complex physical scenario, in which the gas cooling 
time of the forming galaxy is shorter than its dynamical (free-fall) time, so that stars 
form "in flight" , and the subsequent dynamical evolution is a dissipationless collapse. 

The explanation of the observed weak homology of Es is important to understand 
galaxy formation. For example, the presence of a core is usually interpreted as the 
signature of merging of SMBHs, consequence of the merging of the parent galaxies 
(e.g. [332, 333, 243]), while in N-body simulations of repeated equal-mass dissipationless 
mergers the best- fit Sersic index to of the end-products increases with their mass ([236]). 
However, [31, 236] showed that repeated dissipationless merging events fail to reproduce 
the FJ, the Kormendy, and the Mbh-c relations, and also that a substantial number of 
head-on minor mergings make to decrease, bringing the end-products out of the edge-on 
FP (see Sect. 4.1). These results, together with other astrophysical pieces of evidence 
based on stellar population properties such as the Mg 2 -a relation, indicate that dry 
merging cannot have had a major role in the formation of Es, and gaseous dissipation is 
needed (e.g. [253, 260]). 

In alternative (or as a complement) to the merging scenario, it is then of great the- 
oretical interest to explore if (and if so, under what conditions and to what extent) the 
dissipationless collapse of the stellar population produced by a fast episode of gaseous 
dissipation and the consequent burst of star formation is able to reproduce end-products 
with projected density profiles well described by the Sersic law. In particular, following 
the current cosmological picture that galaxies form at peaks of the cold dark matter 
distribution (e.g. [53, 54]), it is natural to investigate dissipationless collapse in two- 
component systems. 

In [331] this process was studied by means of high- resolution two-component N-body 
simulations, in which the collapse of the stellar distribution and the response of the dark 
matter halo are followed in detail. As will be shown, dissipationless collapse in pre- 
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Fig. 28. - Circularized projected stellar density profiles of the end-products of representative one- 
component (left) and two-component (right) collapse events. The dotted lines are the best-fit 
R 1/m models (from [331]). 



existing dark matter haloes is indeed able to reproduce surprisingly well the observed 
weak homology of Es. The flat inner surface brightness profiles of core ellipticals arise 
naturally from dissipationlcss collapse, with the inner core radius i?b determined by the 
coldness of the initial conditions. Two classes of simulations were considered. In the first 
the virialization of a cold, single-component density distribution was followed. In the 
second the initial conditions represent a cold component (stars) deemed to collapse in a 
nearly- virialized live dark matter halo. The initial conditions are the stellar (p*) and the 
halo (ph) density distributions, realized by different combinations of "cold" Plummer and 
7 models (see Sect. 3.2). The corresponding virial ratios /3* and /3h (i.e., the ratio of the 
total kinetic and gravitational energy of the initial conditions), measure the "coldness" 
of the distributions: in a virialized system the virial ratio is 1. 

For a detailed description of how initial conditions are arranged, of the numerical sim- 
ulations, and of the dynamical and structural properties of the end-products see [331]. 
Here it sufficient to recall that the velocity dispersion tensor of the final states is approx- 
imately isotropic in the center and strongly radially anisotropic for r~>ry[ (where tm is 
the spatial half-mass radius), in agreement with previous results (e.g. [60, 330]). 

The projected density profiles of the stellar end-products are fitted over the ra- 
dial range 0.1 <R/R C < 10, which is comparable with or larger than the typical ranges 
spanned by observations ([100]). Apparently all the end-products of one-component col- 
lapse events do not deviate strongly from the i? 1 / 4 law over most of the radial range. 
The end-products of two-component simulations deviate systematically from the i? 1 / 4 , 
and in most cases the profile remains below it at small and large radii. For these systems 
the Sersic index is found to be in the range 1.9<m<12, with average residuals in the 
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Fig. 29. - Circularized projected density profiles of two representative one-component collapse 
simulations with different initial virial ratios. The bars are 1-a uncertainties. The dotted lines 
are the best-fit Sersic models over the radial range 0.1 < R/R c < 10 (from [331]). 



same range as those of one-component collapse events. The quality of the fits is apparent 
in Fig. 28 (left), which plots the surface brightness profile of a projection of one of the 
end-products, together with the best-fit (to = 4.2 ±0.07) Sersic law, and the correspond- 
ing residuals. The average residuals between data and fits are typically 0.04<(A/^i)<0.2, 
where \i — —2.5 logI(R)/(I) e . Figure 28 (right) plots the projected profile of a represen- 
tative two-component simulation with 1Z = 2 together with the best-fit (m = 2.2 ± 0.03) 
model and the residuals. 

So far we have considered the properties of intrinsic and projected density profiles at 
radii >0.1rM (and 0.1i? o )- We now focus on the behaviour of the profiles at smaller radii. 
In fact, the higher resolution of the present simulations allows us to investigate regions 
(down to R <~ 0.01i? o ) comparable to those explored by high-resolution photometry of 
real ellipticals. As apparent from Fig. 29 (top), the end-product density profiles of one- 
component collapse simulations have flat cores at r<0.1rM, in agreement with previous 
studies ([60, 325]). Correspondingly, the projected density profiles are characterized by a 
break radius Rb, in the sense that for R < i?b they stay below the best-fit Sersic profiles 
that matches the profile on the large scale 0.1 < R/R c < 10, as shown in Fig. 29 (bottom 
panels), for a case with /3* = 0.002 (left) and a case with /?* = 0.02 (right). Thus, the 
initial virial ratio /3* determines the radial range over which the final surface density pro- 
file is well fitted by the Sersic law, with colder initial conditions producing smaller cores 
(see Fig. 30). In particular, the fact that the size of the core is correlated (and the max- 
imum central density anti-corrclatcd) with the initial virial ratio is a direct consequence 
of the Liouville Theorem (e.g. [325, 340]). Furthermore, the projected density profiles of 
two-component end-products flatten at small radii and deviate from an inwards extrap- 
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Fig. 30. - Break radius normalized to the effective radius as a function of the initial virial ratio. 
The symbols refer to the average among the three considered projections values, which span the 
range represented by the vertical bars (from [331]). 



olation of the best-fitting Sersic law. However, the flattening is typically more gradual 
and the break is not as apparent as in the one-component cases. This is partly due to 
the fact that the stellar end-products of two-component simulations are characterized, 
in general, by smaller m, corresponding to quite shallow profiles. Note that, according 
to standard interpretations, the central cores observed in several bright ellipticals are a 
consequence of formation through merging, being produced by the interaction of binary 
SMBHs with a stellar cusp (e.g. [243, 332, 333]). However, we have seen that a break 
in the profile at small radii and a flat central core are features produced naturally by 
dissipationless collapse, and this makes dissipationless collapse a plausible alternative to 
the binary SMBHs scenario for the origin of the cores. 

Overall, these results suggest that dissipationless collapse is able to produce systems 
with projected density profiles remarkably similar to the observed surface brightness 
profiles, with high-quality one-parameter Sersic fits even for low m values when non- 
negligible amounts of dark matter are present. In particular, the profiles of the end- 
products of single-component simulations are remarkably similar to those of observed in 
the so-called "core" ellipticals. For instance, the similarity between the profile plotted 
in Fig. 29 (bottom, left panel) with that of the core-elliptical NGC3348, for which the 
reported best-fit is m ~ 3.8 and Rb/R c — 0.016 ([124, 125]), is striking. One can also note 
that, while for m>3 systems the dark matter mass inside the half mass radius is of the 
same order as the visible mass, consistent with observations (e.g. [214, 281], [334] -[339]), 
for very low-m systems dark matter is expected to be dominant. 

5'3. The SLs of dark matter halos. - As anticipated in Sect. 1, much like early-type 
galaxies, also nearby clusters of galaxies define their own FP, a luminosity-radius, and 
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a luminosity- velocity dispersion relations. Understanding similarities and differences of 
SLs in such diverse physical systems can be of the greatest importance in studies of 
galaxy formation, because of the (presumably) different physical processes involved. 

First of all, we recall that well defined SLs are indeed expected for dark matter ha- 
los, on the basis of the simplest model for the formation of structures in an expanding 
Universe, namely the gravitational collapse of density fluctuations in an otherwise ho- 
mogeneous distribution of collisionless dark matter (DM). In fact, the spherical top- hat 
model ([341]) predicts that, at any given epoch, all the existing DM halos have just 
collapsed and virialized, i.e., M = ryCy/G (see Footnote 5). In addition, all the halos 
are characterized by a constant mean density p&, given by the critical density of the 
Universe at that redshift times a factor A depending on z and on the given cosmologi- 
cal model (e.g. [342, 343]). The radius of the sphere containing such a mean density is 
indicated as ta, so that M cx r A . In general, ry ^ r&, but, if for a family of density 
distributions ry/r\ ~ const, then the virial theorem can be rewritten as M cx r\<7 Yl 
and together with M cx r A , it brings to M a Uy, thus providing three relations that 
closely resemble the ones observed for luminous matter. Note that these expectations 
involve the global three-dimensional properties of DM halos, while the quantities entering 
the observed scaling relations are projected on the plane of the sky. However, if DM halos 
are structurally homologous (or weakly homologous) systems and are characterized by 
similar velocity dispersion profiles, as found in cosmological simulations (e.g. [209, 344]- 
[346]), their projected properties are also expected to follow well defined SLs (with some 
scatter due to departures from perfect homology and sphericity). 

Of course, the simple above considerations are not sufficient to explain the observed 
SLs of galaxy clusters, at least for two reasons. The first is that a given potential well 
(as the one associated with the cluster DM distribution) can be filled, in principle, by 
very different distributions of "tracers" (such as the galaxies in the clusters, from which 
the SLs are derived). This means that the very existence of the cluster FP implies a 
remarkable regularity in their formation processes: galaxies must have formed or "fallen" 
in all clusters in a similar way. The second reason is that any trend of the cluster mass- 
to-light ratio (necessary to transform masses, involved in the theoretical relations, into 
luminosities( 7 ), entering the observed ones) must be taken into account for a proper 
interpretation of the observed SLs. 

A distinct but strongly related question about the origin and the meaning of the SLs 
naturally arises when applying the predictions of cosmological models also at galactic 
scales. In fact, while scale-invariant relations are predicted, different slopes of the FJ 
relation are observed for galaxies and for galaxy clusters (see Sect. 2.3). This suggests 
that different processes have been at work in setting or modifying the correlations at 
the two mass scales. As discussed in Sect. 3, the theoretical implications of the scaling 
laws for Es have been extensively explored, and several works have been devoted to 
their study within the framework of the dissipationlcss merging scenario in Newtonian 
dynamics, but surprisingly much less effort has been devoted to the theoretical study 
of the FP of galaxy clusters (e.g. [233]). In order to get a more complete view of the 
problems outlined above, here I report the results obtained in [143], where high-resolution 
N-body simulations have been used to study the scaling relations of very massive DM 
halos. In particular the analysis employed dissipationless simulations with 512 3 particles 



( 7 ) The luminosity of a cluster refers to the sum of the optical luminosities of all its constituent 
galaxies. 
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of 6.86 x 10 10 M Q //i mass each, in a (co-moving) box of side 479 h~ x Mpc. The adopted 
cosmological model is a ACDM Universe with m = 0.3, Oa = 0.7, h = 0.7, spectral 
shape r = 0.21, and normalization to the local cluster abundance, erg = 0.9. From 
this simulation, a sub-sample of 13 halos at z — 0, with masses between 10 14 M Q /h and 
2.3 x 10 15 M Q /h was randomly selected. A first check showed that M o<r\, M«ta4 
with a rms scatter of 0.03 only, and M oc cry 1 , with rms ~ 0.05, for all the halos, as 
expected. 

However, projected quantities are involved in the observations and the first step of 
the analysis is the determination of which (if any) scaling relations are satisfied by the 
DM halos when projected. Therefore the projected radial profiles of the selected halos 
have been constructed by counting the DM particles within concentric shells around the 
center of mass for three arbitrary orthogonal directions, and i?n is defined as the projected 
radius of the circle containing half of the total number of particles. Then, the velocity 
dispersion a^ has been computed from the line-of-sight (barycentric) velocity of all the 
particles within i?n- Since the DM halos (as well as real clusters) are not spherical, such 
a procedure gives different values of i?H and Ch for the three line-of-sights (the maximum 
variations however never exceed 33% and 21% for the two quantities, respectively), so 
that the adopted sample of simulated clusters contains three orthogonal projections for 
each halo. With the projected properties i?H and Oh now available, we have determined 
the best fit relations between M and Rn, and between M and Ch by minimizing the 
distance of the residuals perpendicular to a straight line, and thus obtaining the DM 
analogues of the observed FJ, Kormendy, and FP relations: 

(80) Ma^ 02±al5 , M«4 M14 , 



(81) M OC jR ^±0.05 (7 1.73±0.04. 

in particular, for the last fit it is found that rms = 0.04. Compared to the relations among 
the virial properties, these relations have larger scatter, as expected. The FJ and FP 
have slopes similar to those obtained for the virial quantities, while the M-Ru relation 
appears to be significantly flatter. I stress again that while scaling relations between 
M, r\ (or ry) and ay are expected on theoretical grounds, a tight correlation between 
projected properties is a much less trivial result. In fact, structural and dynamical non- 
homology can, in principle, produce significantly different effective radii and projected 
velocity dispersion profiles for systems characterized by identical M, ry, and ay. It is 
also known that weak homology, coupled with the virial theorem, does indeed produce 
well defined SLs (see Sect. 3). Therefore, the SLs presented in eqs. (80)-(81) are a 
first interesting result. The difference between the values of the exponents appearing in 
eqs. (80)-(81) and those in the virial relations (under the assumption of homology) is the 
direct evidence of weak homology of the halos. Note that this finding is in agreement 
with the results already pointed out by several groups, namely the fact that DM halos 
obtained from numerical N-body simulations in standard cosmologies are characterized 
by significant structural and dynamical weak homology (e.g., [209, 344]). 

Comparison of eqs. (80)-(81) and eqs. (11)-(12) reveals that at the scale of clusters of 
galaxies, the FJ, Kormendy and FP relations of simulated DM halos are characterized 
by different slopes with respect to those derived observationally. What are the impli- 
cations of these differences? In order to answer this question, it is useful to define the 
dimcnsionlcss quantities T = M/L, Q = Ra/R e , and S = a^/a, where R c and a are the 
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quantities related to the optical distribution of galaxies. Focusing first on the edge-on 
FP, from eqs. (12) and (81) one obtains: 

(82) g4m«^°' 42 

Thus, in order to satisfy both the FP of DM halos and the FP of observed clusters 
of galaxies, the product TQ _11 <S -173 must systematically increase as R%' 2 a 0A2 , which 
is approximately proportional to L 3 . In principle, T, Q, and S could all vary in a 
combined and regular way from cluster to cluster, so that eq. (82) is satisfied. Of course, 
given the small scatter around the best fit relation (12), this kind of solution requires 
a remarkable fine tuning of the variations of the three parameters. Alternatively, it is 
possible that only one of the three parameters varies significantly, while the other two 
are approximately constant. This situation is analogous to that faced in the studies of 
the physical origin of the FP tilt of Es, where the so called "orthogonal exploration of 
the parameter space" is often adopted (see Sect. 3). In the present context some of this 
arbitrariness can be removed: in fact, here we assume that 1) the DM distribution in 
real clusters is described by the simulated DM halos, 2) in addition to the edge-on FP, 
we also consider the constraints imposed by the FJ and the Kormcndy relations. These 
two points will allow us to use the orthogonal exploration approach to determine what is 
the most plausible origin of the tilt between the simulated and the observed cluster FP. 

In order to make the DM halos FP reproduce the observed one within the framework 
of the orthogonal exploration approach, we have three different possibilities, each corre- 
sponding to the choice of T, Q, or S as the key parameter, while keeping constant the 
remaining two in the l.h.s. of eq. (82). The two choices based on variations of Q or S 
should be interpreted from an astrophysical point of view as systematic differences in the 
way galaxies populate the cluster DM potential well as a function of the cluster mass. 
However, the orthogonal analysis of the FJ and the Kormendy relations strongly argue 
against these two solutions, since from eqs. (11) and (80) one obtains 

^ ^0.84 £ ^ p0.81 

(83) cx cr , ^Je cx R c . 

Thus, it is apparent that any attempt to reproduce equation (82) by a variation of Q (or 
S) alone will fail at reproducing the FJ (or the Kormendy) relation. In fact, the only 
common parameter appearing in all the eqs. (82) and (83) is the mass-to-light ratio( 8 ) 
T. 

Therefore, while a purely structural (Q) and a purely dynamical (<S) origin of the tilt 
between the DM halos FP and the cluster FP seem to be both ruled out by the above 
arguments, a systematically varying mass-to-light ratio, for Q and S constant, could 
in principle account for all the three considered scaling relations. In particular, from 
eq. (82), T cx L a with a ~ 0.3. Guided by this indication, one can try to superimpose 
the points corresponding to the simulated DM halos to the sample of observed clusters 



( 8 ) Note that the constraints imposed by the FJ and the Kormendy relations should not be 
considered redundant with respect to those imposed by the edge-on FP. These two relations, 
although with a large scatter, describe how galaxies are distributed on the face-on FP. 
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by using T oc M@. It turns out that if 



M 



) 



0.23 



(84) 



T = 280/i 



10 14 M Q //i 



L, 



the edge-on FP of DM halos is practically indistinguishable from that of real clusters (see 
Fig. 5). It is also noticeable (as a non-necessary consequence) that by adopting eq. (84) 
also the face-on FP, the FJ, and the Kormendy relations are very well reproduced, as 
is apparent from Fig. 4. Remarkably, the same trend of the mass-to-light ratio with 
luminosity was found for individual galaxies by the SAURON group (e.g. [281]). It 
would be very interesting to compare the SLs of dark matter halos obtained from high- 
resolution numerical simulations with the SLs of Es and to repeat the above investigation 
in the smaller scale context. 

5'4. Summary. - In this Section three different issues have been addressed. The first 
point concerns the possible simultaneous growth of SMBHs and of the host galaxies at the 
epoch of galaxy formation. The second is about the structural and dynamical properties 
of galaxies formed in dissipationless collapse (the last stages of monolithic-like collapse) 
in pre-existing dark matter halos. The third is the problem posed by the existence of SLs 
of galaxy clusters, which has been discussed in the framework of cosmological simulations 
of structure formation. The main results can be summarized as follows: 

1) At galactic scale, the end-products of one-component simulations of dissipationless 
collapse typically have projected surface brightness profile close to the de Vaucouleurs 
model. When fitted with the Sersic law over the radial range 0.1 <R/R e < 10, the re- 
sulting profiles are characterized by index 3.6<m<8; final states with m>5 are obtained 
only for rather concentrated initial conditions. 

2) The end-products of collapse inside a dark matter halo present significant structural 
non-homology. The best-fit Sersic indices of the stellar projected surface density profile 
span the range 1.9<m<12. Remarkably, the parameter m correlates with the amount of 
dark matter present within R c , being smaller for larger dark-to- visible mass ratios. 

3) The projected stellar density profiles are characterized by a break radius 0.01<i?b/-R e ^0.1 
within which the profile is flatter than the inner extrapolation of the global best-fit Sersic 

law. Colder initial conditions lead to end-products with smaller Rb/R c ; in general, the 
resulting "cores" are better detectable in high- to systems. 

4) For clusters of galaxies, after verifying that DM halos do follow the predictions of 
the spherical collapse model for virialized systems, we have found that also their projected 
properties define a FJ, a Kormendy, and a FP-like relations. However, the slopes of the 
DM halos scaling laws do not coincide with the observed ones, and we have shown that 
the two families of SLs can be reconciled by assuming that the cluster mass-to-light ratio 
T increases as a power law of the luminosity. The required normalization and slope agree 
well with those estimated obscrvationally for real clusters of galaxies. It appears that 
the FJ, Kormendy and FP relations of nearby clusters of galaxies can be explained as the 
result of the cosmological collapse of density fluctuations at the appropriate scales, plus 
a systematic trend of the total mass-to-light ratio with the cluster mass. 

6. Conclusions 



We are finally in the position to connect the different pieces of information described 
in the previous Sections, to see if it is possible to form a plausible scenario in which the 
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existence of the galaxy and central black holes Scaling Laws described in Sect. 2 can be 
traced back to the process of galaxy formation. 

The first important clue is that stellar spheroids are a remarkably regular family 
of stellar systems: the regularity is apparent in terms not only of density profiles, but 
also of orbit composition and of stellar populations. All these indications point towards 
a common formation mechanism, where the galaxy mass has been a major parameter, 
because many galaxy Scaling Laws involve the galaxy luminosity. In Sect. 3 a review 
of the different proposed interpretations of the galaxy Scaling Laws has been presented. 
While a definite answer is not reached yet, it is generally acknowledged that the galaxy 
SLs are mainly due to a systematic variation with the galaxy luminosity of the dark 
matter amount and distribution, of the light distribution (the so-called weak homology) , 
and finally of the metallicity of the bulk of the stellar mass. 

The second piece of information about galaxy Scaling Laws is indirect and comes from 
cosmological simulations: by studying simulations of structure formation on the scale of 
clusters of galaxies, it is found that well defined Scaling Laws are naturally (i.e., by initial 
conditions) imprinted in the resulting dark matter halos. 

Thus, it is tempting to put the two points above together and to suggest that the for- 
mation of stellar spheroids proceeded mainly in a way similar to the monolithic collapse 
scenario. This would explain the galaxy Scaling Laws just as the imprint of the dark 
matter halos Scaling Laws (at the mass scale of galaxies) on the baryons. Numerical sim- 
ulations of fast (dissipationless) collapse in pre-existing dark matter halos can reproduce 
Sersic profiles similar to those observed, from the outer parts of the models down to their 
central regions. Cold dissipationless collapse is a process which is expected to dominate 
the late stages of an initially dissipative process, in which the gas cooling time of the 
forming galaxy is shorter than its dynamical time, so that stars form "in flight" , and the 
subsequent dynamical evolution is dissipationless. Observational evidence supports this 
argument. In fact, the observed color-magnitude and Mg 2 -a relations, and the increase 
of the [a/Fe] ratio with a Q in the stellar population of Es (e.g. [46],??-[349]), suggest 
that star formation in massive ellipticals was not only more efficient than in low-mass 
galaxies, but also that it was faster (i.e., completed before SNIa explosions took place), 
with the time-scales of gas consumption and ejection shorter or comparable to the galaxy 
dynamical time (e.g. [350, 351]) and decreasing with increasing galaxy mass. 

The third piece of information is related to the effects of dry and wet merging on the 
Scaling Laws: in fact we know that ellipticals cannot be originated by parabolic merging 
of low mass spheroids only, even in the presence of substantial gas dissipation (which, 
at variance with dry merging, is able to increase the galaxy central velocity dispersion, 
see also [352])). However, it is also known that SLs such as the FJ, Kormendy, FP, and 
the Mbh-Co relations, when considered over the whole mass range spanned by ellipticals 
in the local Universe, are robust against merging (see also [353]). Thus the galaxy 
Scaling Laws, possibly established at high redshift by the fast collapse in pre-existing 
dark matter halos of gas rich and clumpy stellar distributions (e.g. [354]), can persist 
even in the presence of a (small) number of dry mergers at lower redshift ([355]). If this 
is the case, then monolithic-like collapse at early times and subsequent merging could 
just represent the different phases of galaxy formation (collapse) and evolution (merging, 
in addition to the aging of the stellar populations and related phenomena). 

The possibility that monolithic collapse and successive merging events are just the 
leading physical processes at different times in galaxy evolution, and that they are both 
important for galaxy formation, is perhaps indicated also by a "contradictory" and often 
overlooked peculiarity of massive ellipticals. In fact, while the Kormendy relation dictates 
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that the mean stellar density of galaxies decreases with increasing galaxy mass (a natural 
result of parabolic dry merging), the normalized light profiles of Es becomes steeper and 
their mctallicity increases at increasing galaxy mass (as expected in case of significant 
gas dissipation). Thus, the present-day light profiles of ellipticals could represent the 
fossil evidence of the impact of both processes. 

If the above scenario is correct, then one expects that the star formation history in 
the Universe and the QSO activity should be roughly parallel. It remains to be clarified 
if QSO activity brought the process of galaxy formation to an end, or the star formation 
feedback by ejecting the remaining gas from the galaxy brought the epoch of vigorous 
QSO activity to a conclusion, or finally if a combination of stellar and AGN feedback was 
the key factor. My personal view is that we currently have more indications supporting 
the idea that galaxy formation was stopped more by stellar feedback (i.e., supernova 
heating) than by the AGN feedback (but see [356]). In any case, after the end of the 
fast star formation epoch, necessarily a new evolutionary phase begins for the galaxies 
and their SMBHs. This obvious fact is curiously neglected quite often in the current 
literature, but it is unavoidable. In fact, over a cosmological time in a passively evolving 
galaxy the stellar mass losses amount, over a cosmological time, to a considerable fraction 
(>30%) of the total stellar mass, i.e., ~ 2 orders of magnitude larger than the observed 
SMBHs masses. If only a minor fraction of this recycled gas (the basic ingredient of the 
galaxy "cooling flow" model!) were accreted on the central SMBH, the SLs involving 
the BH masses (such as the Magorrian and the Mbh-Co) would be completely different. 
Thus, the need of an extremely efficient feedback from SMBHs is not required by complex 
physical arguments, but just by the mass budget of the SMBHs. In addition, this feedback 
must be active over the whole galaxy life, and cannot be temporally concentrated just at 
the end of the star formation epoch. However, a moderate accretion from the recycled gas 
by the evolving stars will not destroy an already established SL, such as the Magorrian 
relation, as the available "fuel" for accretion is naturally proportional to the stellar mass 
of the host system. 

I conclude this Review with a brief comment on a recent and very interesting observa- 
tional finding, i.e. the fact that apparently stellar spheroids were much denser than today 
at redshift l<z<2 (e.g. [357, 358, 359], and references therein). The obvious question is 
what mechanism could make a galaxy "expand" . Of course, internal dynamical processes 
cannot be invoked, as their time-scales are measured either by the galaxy dynamical time 
(very short compared to the age of the system) , or by the 2-body relaxation time (which 
is orders of magnitude longer than the age of the Universe). Thus, the only obvious 
possibility is to postulate that a few events of dry merging are common in the life of 
early- type galaxies. This would also help to explain the "central density-slope paradox" 
discussed above. In addition, if dry merging (through a small number of events) is the 
solution to the problem of superdense galaxies then, in order for present-day galaxies 
to obey the Magorrian relation, the SMBHs at the center of the superdense progenitors 
should also follow the same SL, because no significant amount of gas can be accreted on 
the center in a dry merging. Then, the superdense galaxies cannot follow the Mbh-Co 
relation observed in the local Universe because their velocity dispersion is higher than in 
local galaxies of the same mass. In practice, if superdense galaxies are the progenitors of 
the nearby Es, and if their expansion was caused by dry merging, they should obey the 
Magorrian relation observed in the local Universe, but they should fail the local Mbh-c 
relation, by exhibiting a systematically lower zero-point. Testing observationally this 
conjecture would be an interesting goal for the future. 
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